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I. 


Introduction 


In the past several years there has been renewed activity in the 
analysis of the geomagnetic field due to its importance in controlling 
the trajectories of cosmic rays and lower energy charged particles 
trapped in the magnetosphere* Since most of the interest has centered 
on the region from altitudes 100 to 60,000 Km, spherical harmonic repre- 
sentations of the geomagnetic potential have been used almost exclusively. 
Their usefulness in this volume of space is appropriate since the field 
description is smooth and is easy to drop terms from the series as they 
become insignificant with increasing distance from the earth. However, 
there are details in the computation of this field that can confuse the 
novice user of the numerous available computor codes since several different 
conventions are followed. We attempt here to point out these differences / 
in convention and to present a coherent set of programs suitable to dif- 

i 

f erent uses . 1 

II. Formulation 

Making minor changes to the expressions given by Chapman and Bartels 
"Geomagnetism" (1940 p. 639)* we write the potential of the internal field as 

^ n n+1 

V = a H (r ) (g n,m cos m tP h * sin m<P ) p n » m (q) 

n=l m=0 I I 

where a - radius of the earth 

r = geocentric distance 
g,h=Gauss coefficients 
^ = longitude 

* We will hereafter refer to this book as GM 
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© = colatitude 

p n,m (q) = associated Legendre functions (Gauss normalized) 

from which three orthogonal components may be derived by taking the gradient 
B = + 7 V to give 

n n+ 2 ti m n m 

Bq = AJL =2 2 /a\ (g n,m cos myH-h * sin mj 0 m) j£_P ’ 

r © n=l m=0 ' r / d 9 

B = l c>V=2 2 /av l+2 -JL (~g n,m sin + B n,m cos mf) p n,m (6) 

¥ r sin © b f n=l m=0 \r) sin© 

B = _^_V = ^ /af +2 (n+1) (g n,m cos mf+h n,m sin mf )P n ’ m (©) 

r dr n=l m=6 ' r' 

i-5 2 2 

and B = y Bq + B ^ + B r 

We here follow the notation in GM pp. 610-611 in denoting the (Gauss- 
Laplace) funtions P n,m (©) which are distinct in normalization from 
the Neumann functions P n m (©) an B the Schmidt functions P n (©) . The re- 
lations between these functions is given explicitly in GM pp. 610—611, 

Equations 11, 19 said 20. Here it is useful to note that 

(2n)!/2 n n! = (2n-l) ! ! 

The generating functions used to compute P and _3_P are as 

3 © 

follows: 

p°>° = i yp°’° = ° 

? © 

P n > n = (sin©) P 11 " 1 ’ n - X 

d_P n > n = (sin©) 2l n ~ 1,n ~ 1 + (cos©) P n_1 ’ n_1 ( n 9 1} 

d 0 3 6 

and for m 4 n 2; 1 
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p n > m = (cos©) pn-l,m _ K n,m p n-2,m 

d P n ,m = (cos e) ^ pn -l,m _(s jn Q) pn-l,m \ pn-2,m where 

be ie jg 

K n,m = (n-l) 2 -m 2 n >l 

(2n-l)(2n-3) 

K n,m =0. n = 1 

For the dipole term of degree n=l and order m=o the previous 
expressions can be used to derive the field at the earths surface 
(r=a) as follows: 

Bq = g 1 ’ 0 ll 1 ’ 0 = - g 1 ’ 0 S inQ 

£ e 

B ^ =0 

B r = -2 g 1,D P 1 ’ 0 = -2g 1>0 Cos e 

Thus for g 1 * 0 ^>0, Bq < 0 for all 9 and B p < 0 in the northern hemisphere 
(e •<90°). These directions are different from those conventionally 
used in geomagnetism as illustrated in Figure 1 (taken from GM pg. 2). 

In this figure the three orthogonal components of the field are X (=n'orth) , 
Y (= east) and g (= downward) . Except for a slight difference of 
definition between geodetic and geocentric coordinates to be discussed 
subsequently, Bq = -X, B r = -g and B^ = Y. We realize that in choosing 
to perform the evaluation of the field components utilizing P n > m (9), 
instead of the Schmidt quasi-normalized functions P^ (6) , we are creating 
a problem, in that by international convention the g's and h's for the 
expression of the geomagnetic field are to be used with the PjJ's. Comments 
on this point are to be found in GM pg. 638 as follows: 
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"Workers in this field are advised to use the partly normalized 
functions p£} introduced by Adolf Schmidt. They have the ad- 

vantage that the mean square value of P™ is the same as that for 
P®, so that the coefficients in the series indicate the relative 
importance of the various terms. This choice would also remedy, 
for the future, the regrettable confusion of notation which has 
existed in geomagnetic literature in the past. A resolution re- 
commending the use of Schmidt's functions was adopted by the 
International Association of Terrestrial Magnetism and Electricity 
at its Washington meeting in 1939" • 

The reason that the gauss normalization is adopted for computer 
use is that in the evaluation of Bq, and Bj. the choice of unity for 
the first term of the generating function saves two multiplications 
for each n and m. Since computer codes are being used extensively, this 
slight shortening of the computation time is considered worth the 
complication of converting g^, h™ to g n,m , h n,m . We maintain the 
difference of convention used by Finch and Leaton (I960) whereby there 
is a change of sign between the two types of coefficients and a corre- 
sponding difference either in the sign of the potential V or in taking 
the gradient V V to compute the field. Thus for the earth's field the 
dipole term > 0 whereas gq <-0. 

The factors Sa m used for converting the "Schmidt normalized" co- 
efficients g^, h^, to the "Gauss normalized" g 11 ’ 111 , h n,m are as follows: 



S n,m = S n,m-1 J- 

T n + m 


where J = 2 for m = 1 
J = 1 for m > 1 
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Oblate Earth 

The above formulation is rigorously correct only for a sphere. 

In the past geomagneticians have chosen to ignore the fact that the 
earth is non-spherical even to the extent of making no allowance for 
its oblateness. To again quote GM (pg. 641 ): 

"While it is of considerable theoretical interest to consider 
the influence of the earth's oblate form, it must be admitted that 
this refinement has not added much to our actual knowledge of the 
field. It makes no serious difference to the estimates of either 
the exterior part, the interior part, or the non-potential part 
of the results. No significant physical property of the geo- 
magnetic field can therefore be attributed to the earth's oblateness." 
However, as the accuracy of evaluation of the earth's field increases 
it is obvious that it will eventually be essential to take the earth's 
true shape into account. So long as the evaluation of the harmonic 
coefficients is done in spherical coordinates, r,9 and 1 ^, the resulting 
fields B r , Bq and will be in strict geocentric directions. The 
only constant pertaining to the earth is the radius a used in the potential 
function. To be consistent with past determinations we have adopted the 
mean value of 6371.2 Km for a. With this philosophy in mind the only 
problem is that of converting positions in geodetic coordinates to 
geocentric before evaluating the field and of transforming the resulting 
geocentric field vectors back to the geodetic system to compare with 
data measured in this system. The expressions for this conversion are 
given here only for the oblateness of the earth. Referring to Figure 2 
we may write : 
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tan'X 


and R^ = +• 



cos^X +■ 


B* 


sin^A + 


A^~ cos^ X + B^ sin^ A 

A2 COS^ A + b 2 sin2 X 


where : 

h = height above the geoid 
= geocentric latitude (90°-9) 

A = geodetic latitude 
R = geocentric distance 
A = equatorial radius (=6378.165 Km) 

B = Polar radius (=6356.783 Km) 

Using A and h the geocentric quantities 9= 90 - ip and R may thus 
be calculated for use in the previous expressions for B. The conversion 
from Bq, B r to X and Y can then be done by the rotation: 
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X = -Bq cos S -Bj. sinS 
g = Bq sin<£ -B r cos £ 
where = TV - 

For the use oi coefficients of the field derived from surface 

data using the approximation that the earth is spherical, the proper 

\ 

relations are: 


6 = 90 - A 

X = ~B 0 

r = 6371.2 + h. 

g = -B r 


III. General Purpose Field Program 

The purpose of the first set of programs is the evaluation of the 
geomagnetic field components X, Y, 3, F using any available set of 
harmonic coefficients at any point in normal geodetic coordinates 
Latitude, Longitude and Altitude. 

The logic of the main control program FIELDG is given in 
Figure 3 and the FORTRAN listing in Table I. As can be seen in the 
listing, the variables are: 

INPUT: DLAT = Geodetic latitude Ain degrees (north positive). 

DLONG = East longitude in degrees. 

ALT = Height above the earth in Km. 

TM = Time in years for which field is desired (e.g. July 

15, 1964 = 1964.54) 

NMX = Maximum degree N (=n + l) of spherical harmonic 
coefficients to be used. 


L = Code to read in new coefficients when > 0. 
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OUTPUT ; X = North geodetic component (=BN)* 

Y = East component (=BE)* 
g = Vertical component (=BV)* 

F = Total field (=B)* 

The units of the output field components will be the same as those 

of the coefficients. 

The logical blocks of the program are thus: 

(1) Set Up Constants - 

On the first call of the subroutine the constants are 
set up for later use. The earth ellipsoid used here is 
that discussed by Kaula (1963), with an equatorial radius 
of 6378.165 Km and flattening factor of 1/298.3. (Polar 
radius = 6356.783 Km) 

(2) Coefficient Set up - 
(a.) L Test 

This block branching from L < 0 reads a set of coefficients 
from cards in the format given in section V of this report. 

The L value is then internally set 0 if time terms are pre- 
sent or negative if not, for later use by the subroutine. 

This test is made only on GT (N,M) on the assumption that 
at least one such time derivative will be present. L should 


*The labels BN, BE, BV, B were previously used in the Vanguard 3 
Report ( Cain et al, 1962) and some subsequent programs. 
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only be set positive by the user if a new set of coeffi- 
cients is to be read in. Since L is an input parameter 
altered by the subroutine it should always be given as a 
variable in the calling sequence rather than a fixed point 
constant. 

(b.) K Test 

The branch on K=0 is the conversion of the Schmidt 
to Gauss normalized coefficients. For K^O the coefficients 
are assumed to be Gauss normalized and not converted. 

(c.) Compute Coefficients for a new Time 

The coefficients TG and TH are computed from the 
G's and H's after being read in (L > 0) and thereafter 
if L=0 and the time variable (TM) changes. Since the 
secular change is small, many users may only wish to 
change TM in increments of the order of a year. 

(3) Computations 

The computations from 'Comp 1' up to where the FIELD sub- 
routine is called calculate the proper inputs for FIELD. The 
NMAX used in FIELD is limited to the smallest of MAXN or NMX. 

(4) Field Subroutine 

This subroutine is given here in two versions as listed in 
Tables II and IV respectively. Although the second version is 
discussed under the next heading, it is so constructed that it 
can also be called by F IEL DG. The version in Table II is only a 
slight modification of that given in the Vanguard 3 report ( Cain 

et al 1962). The variables are: 
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INPUT: ST = sine, CT = cosO 

SPH = sin y>, CPH = cos y> 

R = geocentric distance (Km) 

NMAX - (n + l) maximum 
OUTPUT: BT = B 0 , BP = BR = Bp, B 

The gauss normalized coefficients G and H may be fed from the TG 
and TH in FIELDG through COMMON. In the event that one wished 
to use Schmidt normalized coefficients directly instead of con- 
verting it would only be necessary to add the factor SHMIT (N,M) 
from the FIELDG program (after proper set-up) to the statements 
PNM = P (N,M) * AR and TEMP = G (N,M) * CP(M) + H(N,M) * SP (M) . 
It is these NMAX -1 pairs of multiplications that are saved on 
each subroutine call by the use of gauss normalized coefficients. 

One should note that although FIELD is here given as a sub- 
program of FIELDG, it can also appropriately be used independently 
provided the input variables and coefficients are in the proper 
form. 

IV * Special Purpose Field Programs 

(a) FIELD (non-indexed version): 

For ultimate economy a machine-language hand coded subroutine 
is the most efficient. However, since the use of subscripted 
variables also increases execution time, particularly on some com- 
puter compilers which do not handle indexed variables very effi- 
ciently, a FIELD program with constants or constant indices could 


shorten the execution time. This second version of FIELD is 
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provided with the intent of supplying this extra efficiency. 
Since it is a fairly lengthy program, it is expedient to 
generate the FORTRAN text by the program code given in Table 
III. This FIELD generator is set to generate the text to 
NMAX = 18 but of course can be cut down to meet any specific 
NMAX with a resulting savings in core locations. The text 
actually listed in Table IV is a somewhat more compressed 
version of that which the program generator produces. How- 
ever, the compiled version of the two programs are identical, 
(b) Timing: 

An approximate timing for the FIELDG subroutine for NMAX=8 
and TM changing is 5*5 ms on an IBM 7094 Model I. If TM 
remains constant the execution time drops to about 2.5 ms. 

A table of the execution time for the two versions of FIELD 
on various machines is as follows (NMAX = 8) : 

FIELD TIME (millisec) 



indexed 

non- indexed 

IBM 7094 (II) 

6.8 

4.4 

CDC 3600 

8.7 

6. 

IBM 7094 (I) 

12.4 

8.3 

UNIVAC 1107 

16 . 9 

12.6 


IBM 1410 


3700. 


2500. 
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(c) Decrease of NMAX with h 

Although the dropping of terms from the expansion with 
increasing r must be viewed according to the special use of 
the program, we include here two FORTRAN statements that may 
be included to cut execution time for large R: 

NMAX = MAXI (3., ABS (7.5/ALOG (1.001 + ALT/6371.2)* 
NMAX = MINO (NMAX, MAXN) 

where MAXN = maximum number of coefficients available. 

Instead of performing the above computation each time 
one could instead of course refer to an approximate table of 
values of R for FIELD or ALT (Km) for FIELDG as follows: 

NMAX ALT R/6371.2 

3 71000 12.0 

4 35000 6.5 

5 22000 4.5 

6 16000 3*5 

7 12000 2.9 

8 10000 2.6 

9 8000 2.3 

10 7000 2.1 

(d) Coefficient converter 

For comparison of coefficients of different normalization 
it may be useful to have a simple subroutine for coefficient 
conversion. Such a subroutine is given in Table V with for- 
mulation similar to that used in FIELDG. 


^Fortran IV AL0G is log e 
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V. Spherical Harmonic Coefficients 

Collected herein are some of the published sets of spherical 
harmonic coefficients in the format written by FIELDG in units of 
gamma (=10 _ ^Gauss) . (Table VI) 

VI. Sample Values 

A few sample values of field generated by the FIELDG program 
using the set of parameters marked "April 64 coefficients" are as 
follows: 

TM = 1960.0 (=Epoch) , NMX = 8 


DLAT 

DLONG 

ALT 

X 

Y 

2 

F 

-60 

-180 

0 

10362 

8221 

-64160 

65509 

-60 

-180 

100 

9822 

7826 

-61048 

62327 

-60 

0 

0 

16356 

-6532 

-31525 

36111 

0 

0 

0 

28047 

-5969 

-11352 

30840 

0 

0 

100 

26710 

-5695 

-10446 

29240 

30 

-60 

0 

22650 

-6341 

43057 

49063 

60 

120 

0 

14060 

-2976 

59558 

61268 
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SUBROUTINE FIELDG 



^ Any 
Time Ter ms 7 s * 

13 

Yes 

L = 

0 

12 < 

s' 

> No 




L « 

-L 




14, 





\ 





=0 

(Schmidt) 


DLAT = Degrees Latitude 
DLONG - Degrees East Longitude 
ALT = Altitude (Kilometers) 

TM - Time in years 
NMX = Maximum degree + 1 

L Control variable (Set > 0 to read in 

coefficients) 

X = North Component of Field 
Y - East Component 

Z - Vertical Component (Positive Down) 

F = Total Field 

VARIABLES DESCRIBING INPUT COEFFICIENTS: 

J f 0 Coefficients derived using spherical 
earth approximation 

J * 0 Coefficients derived using oblate 
earth 

K = 0 Schmidt Normalized 
K / 0 Gauss Normalized 


♦Numbers at Branch Points are statement numbers 


Figure J 



Table I 


SUBROUTINE F I ELDG ( DL AT , DL O.MG, ALT , TM, NMX , L , X , Y , 7 , F ) l 

c geomagnetic field using any set of coefficients 2 

C L POSITIVE READS A NEW SET OF COEFFICIENTS. FIRST CARD 3 

C READ AS IDEM IFICAT ION IN COLUMNS 10-72. IF COL. 1 ZERO THEN A 

C COEFFICIENTS ARE FOR OBLATE EARTH, IF NON ZERO THEN SPHERICAL 5 

C EARTH. IF COL. 2 ZERO THEN SCHMIDT NORMALIZED, IF NON- ZERO THEN 6 

C GAUSS NORMALIZED. COLS. 4-9 CONTAIN EPOCH YEAR. 7 

C PLAT, ALT, X AND Z ARE GEODETIC 8 

C REOU IRES SUBROUTINE F I ELD{ ST , CT , SPH , CPH , R , NM AX, R T , BP ,BR , B ) 9 

DIMENSION Hi 16,18), GI 18,18) »HT( 18,I8),GT(18,18),HTTI18,I8),GTT(18, 10 

U3),SHM1T( 18, 18), AID! ID II 

COMMON /C0EFFS/TGI18, 18),TH(13,13) 12 

IF (SHM1T{1,1)+1. ) 1,3,1 13 

A=6378. 165 14 

FLAT=1.— 1./298. 3 15 

A2=A**2 16 

A4=A**4 17 

32= ( A»FLAT )**2 18 

A2B2=A2*( l.-FLAT**2) 19 

A4B4=A4* ( l.-FLAT**4) 20 

SHM I T ( 1 , 1 ) =— 1 . 21 

DO 2 N=2 » 1 8 22 

SHM IT (N,1)=SHMIT(N-1, I )*FLOAT { 2*N-3 ) /FEOATIN-1 ) 23 

J = 2 24 

DO 2 M=2,N 25 

SHMIT(N,M)=SHMIT(N,M-1 )*SQRT( FLOAT! ( N-M+ 1 )* J ) /F LOAT { N+M-2 I) 26 

2 J=i 27 

3 IF (L) 19,4,5 28 

4 IF (TM-TLAST) 17,19,17 29 

5 READ (5,6) J , K, TZERO, ( A I0( I ) , 1=1, 1 1 ) 30 

6 FORMAT (211, IX , F 6 . 1 , 10A6 , A 3 ) 31 

WRITE (6,7) J,K,TZEROt ( AID! I), 1=1,11) 32 

7 FORMAT (2I3,5X,6HEP0CH=,F7.1,5X10A6,A3) 33 

MAXN=0 34 

TEMP=0. 35 

8 READ (5,9) N, M, GNM, HNM , GTNM, HTNM , GTTNM, HTTNM 36 

9 FORMAT (213, 6F 11. 4) 37 

IF (N) 11,11,10 38 

10 M A X N = ( M A X 0 ( N » M AX N ) ) 39 

G(N,M)=GNM 40 

H ( N, M ) =HNM 41 

GT ( N, M ) =GTNM 42 

HT (N,M)-=HT NM 43 

GTT(N,M)=GTTNM 44 

HTT ( N,M) =HTTNM 45 

TEMP=AMAX1( TEMP,ABS(GTNM) ) 46 

GO TO 8 47 

11 WRITE (6,25) ((N,M,G(N,M).H(N,M),GT(N,M),HT(N,M),GTT(N,M) ,HTT(N,M) 48 

I,M=1,N) ,N=2,MAXN) 49 

IF (TEMP) 13,12,13 50 

12 L=— L 51 

GO TO 14 52 

13 L=0 53 

14 IF (K) 17,15,17 54 

13 DO 16 N=2 , MAXN 55 

DU 16 M=1 , N 56 

G(N,M)=G(N,M)*SHMIT(N,M) 57 

H ( N , M ) =H ( N, M ) * SHM I T ( N , M ) 58 

GT ( N , M ) =GT ( N , M ) *SHM IT ( N , M ) 59 

HT(N,M)=HT(N, M)*SHMIT(N,M) 60 

GT T ( N » M ) =G T T ( N » M ) * SHM I T ( N , M ) 61 


Table I Cont. 

16 HTT(N,M)=HT1 {N,M)*SHM IT(N,M) 62 

17 r = TM-TZ ERU 6 3 

00 13 M = 2,MAX\I 64 

00 IB M=l»N 65 

1 G ( .N » M ) -G( N , M ) + T* ( G T ( N , M ) -i-G I T ( N , M ) * T ) 66 

13 TH(M,M)=H(N,M)+T*(HT(.'B,M)+HTT(N,H)*T) 67 

T L AS r = T M 68 

19 SINLA=S INI CL AT/ 57. 295 7755) 69 

RLUNG = r>L0NG/57. 2957795 70 

CRH=CUS I RLONG ) 71 

$l>H=S INIRLONG) 72 

IF (J) 20,21,20 73 

20 R = ALT + 6371 .2 74 

C T = S I N L A 75 

GO TO 22 76 

21 S I NL A2 = $ I NL A* *2 77 

COSLA2=l.-SINLA2 78 

DEN?=A2-A2B2*SINLA2 79 

DEN=SQRT (DEN? ) 80 

F A C = ( ( (ALT*UEN)+A?)/( ( ALTOEN) +B2 ) ) **2 81 

CT-SINLA/SQRI l F A C * C 0 S L A 2 + S IN L A 2 ) 82 

R=SQRT ( AL1 *( ALT + 2.*DEil)+(A4-A4B4*SINLA2)/DEN2) 83 

22 S T = S 9RT ( 1 . — C T * *2 ) 84 

N!-1AX= M 1 00 ( NMX , M A X N ) 86 

CALL FIELD ( S T ,CT , SPH , CRH , R , MM AX , B T , Y , BR , F ) 86 

IF (J) 23,2 4,23 87 

23 X = -BT 88 

7 = -BR 89 

RETURN 90 

C TRANSFORMS FIELD TO GEODETIC DIRECTIONS 91 

24 S I NO = S INL A*ST— SORT ( COSL A2 } »CT 92 

C OSD= SQRT ( 1.0— SIND* »2) 93 

X = -BT *CGSD— BR * SIND 94 

Z=BI *SIND— BR*COSD 95 

25 I DRIB AT (6H0 N M6X 1HG10X 1HH9X2HG T9X 2HHT 9X 3HG I T8X3HHT T // ( 2 I 3 , 6F 1 1 . 4 96 

1 ) ) 97 

RETURN 98 

LNC 99 



Table II 

SUBROUTINE FIELD ( S l , CT , SPH ,CPH , R , NMAX, BT, BP , BR, B ) 

COMMON /CG£FFS/G( 18, 18 ) ,H< 13,18) 

01 MENS ION P( 18, 18 ), DP l 18, 1 8), CONST ( 18 , 1 8 ) , SP ( 18) ,CP <181 ,FN( 18) , FM { 
118) 

IF <P(1,1)-1.0> 1,3,1 

1 P(l,l)=l. 

DPI l,l)=0. 

5P ( 1 ) =0 • 

CP ( 1 ) =1 • 

DO 2 N=2 , 1 6 
FN(N)-N 
DO 2 M= 1 , N 
FM ( M ) =M— 1 

2 CONST (N,M)=FLOAT ( ( N— 2 ) **2~I M— 1 ) ** 2 ) /FLOAT { 1 2*N-3 ) * { 2*N-5 ) ) 

CONS T ( 2 ♦ 2 )=0 • 

3 SP ( 2 ) = SPH 
CP ( 2 ) = CPH 

00 4 M=3 , NMAX 

SP (M)=SPI 2)*CP(M-l)+CP(2)*SP(M-l ) 

4 CP{ M)=CP(2)*CP{M-1 )-SP(2)*SP(M-l) 

A0R = 6 37 I * 2 /R 

AR-A0R*«2 
B1 =0. 

BP = C. 

BR =0* 

00 8 N=2 , NMAX 
AR- AQR*AR 
DO 8 M- 1 ,N 
IF CN-M) 6,5,6 

5 P(N,N)=ST*P(N-I,N-1) 

DP(N,N)=ST*DP(N-1,N-1)+CT*P(N-1,N-1 ) 

GO TO 7 

6 P(N, M)=CT*P(N-1,M)-C0NST(N,M)*P(N-2,M) 

DP(N»M)=CT*DP{ N— I , M ) —ST *P ( N- L , M)— CONST { N , M ) *0P ( N-2 , M ) 

7 PAR= P ( N , M ) * AR 

TEMP- G (N, M ) *CP ( M ) +H ( M , M ) *SP ( M ) 

8T=BT+TEMP*DP(N,M)*AR 

BP- BP— ( G { N , M ) *SP( M) — H ( N, M ) *CP ( M ) )*FM{M)*PAR 

8 BR=BR-TEMP*FNI N)*PAR 
BP=3P/ST 

B=SQRT(BT*BT+BP*BP+BR*8R) 

RETURN 
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Table III 


C PUNCHES FORTRAN TEXT FOR NON-SUBSCRI PTEO FIELD PROGRAM 

PUNCH l 4 

1 FORMAT (6X.50HSUBR0UTINE FIELD I ST , CT , SPH , C PH , R, NMAX, BT ,BP , BR, B ) ) 5 

PUNCH 2 6 

2 FORMAT (6X,3IHCOMMON/COEFFS/G( 18,18), H(18, 18)) 7 

PUNCH 3 8 

3 FORMAT ( 6X,11HAR=6371.2/R) 9 

PUNCH 4 10 

4 FORMAT { 1HC,59X,4HN= 2) 11 

PUNCH 5 12 

5 FORMAT ( 6X,7HSP2=SPH ) 13 

PUNCH 6 14 

6 FORMAT (6X,7HCP2=CPH) 15 

PUNCH 7 16 

7 FORMAT ( 6X,6HP21=CT ) 17 

PUNCH 8 18 

8 FORMAT ( 6X,6HP22= ST ) 19 

PUNCH 9 20 

9 FORMAT (6X,9HDP21=— P22) 21 

PUNCH 10 22 

10 FORMAT <6X,8H0P22=P21 ) 23 

PUNCH 11 24 

11 FORMAT { 6X, 12HA0R=AR*AR*AR ) 25 

PUNCH 12 26 

12 FORMAT (6X,24HC2=G(2,2)*CP2*H(2,2)*SP2) 27 

PUNCH 13 28 

13 FORMAT (6X,33HBR=-(A0R+A0R)*(G(2,1)*P21+C2*P22) ) 29 

PUNCH 14 30 

14 FORMAT ( 6X , 2 9HBT =AOR*lG( 2, 1 ) *DP2 l+C2*DP22 ) ) 31 

PUNCH 15 32 

15 FORMAT (6X.35HBP = AOR* { H( 2,2) *CP2-G( 2,2)«SP2)*P22) 33 

00 57 N=3, 1C 34 

NM=N-1 35 

NM2=N-2 36 

PUNCH 16,NM»N 37 

16 FORMAT (6X,8HIF(NMAX-, 12, 6H) 1,1, ,12) 38 

PUNCH 17, N 39 

17 FORMAT ( 1HC, 59X,2HN=, 12 ) 40 

IF t MODI N, 2) ) 18,21,18 41 

18 N02=NM/2+i 42 

PUNCH 19,N,N»N02,N02» N02 43 

19 FORMAT ( I2,4X,2HSP,I2,4H=(SP,I2,3H+SP,I2,4H)*CP,I2) 44 

PUNCH 20,N,M02»N02, N02 , N02 45 

20 FORMAT ( 6X , 2HC P, I 2 , 4H= ( CP , 1 2 , 3H+SP , 1 2 , 5H ) * ( CP , I 2 , 3H-SP , 12 , 1H ) ) 46 

GO TO 24 47 

21 PUNCH 22,N,N,NM,NM 48 

22 FORMAT ( 12 ,4X, 2HSP , I 2 , 7H=SP2*CP, I 2 , 7H+CP2*SP , 1 2 ) 49 

PUNCH 23,N,NM, NM 50 

23 FORMAT (6X,2HCP,I2,7H=CP2*CP,I2,7H-SP2*SP,I2) 51 

24 IF (N-3) 32,25,32 52 

25 PUNCH 26 53 

26 FORMAT ( 6X, 23HP3 l=P2 1*P 2 1-0 . 333333333 ) 54 

PUNCH 27 55 

27 FORMAT ( 6X , 1 1HP32=P2 1*P22 ) 56 

PUNCH 28 57 

28 FORMAT ( 6X , 1 1HP33=P22*P22 ) 58 

PUNCH 29 59 

29 FORMAT ( 6X , 13HDP31=-P32-P32 ) 60 

PUNCH 30 61 

30 FORMAT ( 6X , 16HDP32=P2 1*P2 1-P33 ) 62 

PUNCH 31 63 
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31 FORMAT l 6X, I OH DP 3 3= -OP 31 ) 64 

GO TO 43 65 

32 DO 40 M=1,NM 66 

I NUM= ( N+M-3 ) * ( N— M— 1 ) 67 

I DEN= ( N+N— 3 ) * ( N+N— 5) 68 

1 C 1= I 10000*INUM)/IDEN+10000 69 

INUM=MOD( 10000*INUM, I DEN) 70 

IC2=(10000«I NUM )/IDEN +10000 71 

PUNCH 33,N,M,NM,M 72 

33 FORMAT (6X, 1HP , 2 1 2, 6H=P21*P » 2 1 2 ) 73 

IF (M-NM) 34,36,34 74 

34 PUNCH 35, IC1, IC2.NM2.M 75 

35 FORMAT ( 5X ,4HS-0 . , 2 14 , 2H*P , 2 12 ) 76 

36 PUNCH 37»N,M,NM,M,NM,M 77 

37 FORMAT ( 6X , 2 HDP , 2 1 2 , 7H=P2 1 *DP , 2 12 , 7H+DP21 *P , 21 2 ) 78 

IF (M-NM) 38,40,38 79 

38 PUNCH 39, IC1, IC2,NM2*M 80 

39 FORMAT ( 5X , 4H$-0 . , 2 1 4 , 3H*OP , 2 12 ) 81 

40 CONTINUE 82 

PUNCH 4 1 , N , N » NM , NM 83 

41 FORMAT (6X»1HP,2I2,6H=P22*P»2I2) 84 

PUNCH 42 »N,N,NM, N,NM 85 

42 FORMAT ( 6X , 2H0P , 21 2 , 1H= , I 2 , 4H. 0*P , 2 12 ) 86 

43 PUNCH 44 87 

44 FORMAT ( 6X , 1 0HA0R=A0R» AR ) 88 

00 46 M=2,N- 89 

PUNCH 45,M,N,M,M,N«M,M 90 

45 FORMAT ( 6X, 1HC , 12 , 3H=G ( , I 2 , IH, , 12, 4H ) »CP, 12 , 3H+H< , I 2 , 1H , , 12 , 4H ) *SP 91 

1,12) 92 

46 CONTINUE 93 

PUNCH 47 , N, N , N 94 

47 FORMAT (6X,8HBR =BR - , 1 2 , 10H. 0*A0R* ( G ( , 1 2 , 5H, l ) *P » 1 2 » 1H1 ) 95 

DO 49 M=2 , N 96 

PUNCH 48,M,N,M 97 

48 FORMAT (5X,3H$+C, I2»2H*P»2I2) 98 

49 CONTINUE 99 

PUNCH 50 100 

50 FORMAT (5X,2HS)) 101 

PUNCH 51 ,N,N 102 

51 FORMAT ( 6X , 1 5HBT =BT +AOR* ( G(, I2,6H,1)*0P, I2.1HI) 103 

00 53 M=2,N 104 

PUNCH 52 , M,N, M 105 

52 FORMAT ( 5X,3H$+C, 12, 3H*0P,2 12) 106 

53 CONTINUE 107 

PUNCH 50 108 

PUNCH 54 « N,N,N 109 

54 FORMAT (6X.16HBP *BP -AOR* < ( G( , I 2 , 10H , 2 ) *SP2-H ( , I 2 , 10H , 2 ) *CP2 ) *P , I 110 

12 , 1H2 ) 111 

DO 56 M=3, N 112 

MM=M— l 113 

PUNCH 55,MM,N,M,M,N,M,M,N,M 114 

55 FORMAT (5X,2Hi+,I2,6H.0*(G(,I2,lH,,I2,4H)*SP,I2,3H-H(,I2,lH,,I2,4H 115 

1)*CP, I2,3H)*P,2I2) 116 

56 CONTINUE 117 

PUNCH 50 118 

57 CONTINUE 119 

PUNCH 58 120 

58 FORMAT ( 1H1 , 5X9HBP=BP/P2 2) 121 

PUNCH 59 122 

59 FORMAT (6X,25HB=SQRT( BR*BR+BT*BT+BP*BP) ) 123 

PUNCH 60 124 




Table IV 


SUBROUTINE F IFLD ( ST , CT , SPH, CRH, R , NMAX , B T , BP , BR , B ) 1 

COMMON/COE FFS/GI 18,18),H( 18, 13) 2 

AR=637i.2/R 3 

C N =24 

1 SP2=SPH 5 

CP2=CPH 6 

P21=CT 7 

P22=ST 8 

OP2 I=-P22 9 

DP22=P2I 10 

AOR=AR*AR*AR 11 

C2=G( 2, 2 ) *CP2+H( 2, 2 )*SP2 12 

BR=-( AOR+AOR )*(G(2»1)*P21+C2*P22) 13 

BT=AOR*(G<2,l)*DP21+C2»DP22) 14 

BP=AUR«(HI2»2) *CP2— G( 2, 2) *SP2 ) *P22 15 

IF ( NMAX- 2) 18,18,2 16 

C N= 3 17 

2 SP3=(SP2+SP2)*CP2 18 

CP3=(CP2*SP2)*( CP2-SP2 ) 19 

P3l=P2l*P21-0. 333333333 20 

P32=P21*P22 21 

P33=P22*P22 22 

0P31=— P32— P32 23 

DP32=P2l*P21-P33 24 

0P33=— 0P31 25 

AOR=AOR*AR 26 

C2=G(3,2)*CP2+HI3,2)*SP2 27 

C3=G(3,3)*CP3*H(3,3)*SP3 28 

9R=3R-3.0»A0R»IG( 3, 1 )*P31+C2»P32+C3*P33) 29 

B l=BT ♦ AOR» (G( 3,1)*0P31+C2*DP32+C3»0P33) 30 

BP=BP-AOR*( (G< 3,2 )*SP2-H( 3,2 )*CP2)*P 32 +2.0*1 G( 3,3) *SP3-H( 3,3) *CP3) 31 

l*P33) 32 

IF ( NMAX— 3) 18,18,3 33 

C N= 4 34 

3 SP4=SP2»CP3+CP2*SP3 35 

CP4=CP2*CP3-SP2»SP3 36 

P4 1 = P2 1 *P3 1— 0 . 2 66 6666 6* P 2 1 37 

DP4l=P21*DP3l+DP21*P3l-0.26666666*DP21 38 

P42=P2 1 *P32-0. 20000000*P22 39 

DP42=P21*DP32+DP21*P32-0. 20 000000* DP 2 2 40 

P43=P21*P33 41 

DP43=P21*DP33+DP21*P33 42 

P44=P22*P33 43 

DP44=3.0*P43 44 

AOR=AOR*AR 45 

C2=G(4,2)*CP2+H(4,2)*SP2 46 

C3=G(4,3)*CP3+H(4,3)*SP3 47 

C4=G(4,4)*CP4+H(4,4)*SP4 48 

8R=BR-4.0*AQR*(G( 4, l) *P41+C2*P42+C3*P43+C4*P44) 49 

8T=BT+A0R*(G(4,1 )*0P41+C2*DP42+C3*DP43+C4*DP44) 50 

BP=BP-AUR*t (G(4,2)*SP2-H(4,2)»CP2)*P42+2.0*(G(4,3)*SP3-H(4,3)*CP3) 51 

1* P43+3.0 *( G ( 4, 4 ) * SP4-HI 4 , 4 )*CP4 ) *P44) 52 

I F ( NMAX— 4) 18,18,4 53 

C N= 5 54 

4 SP5=l SP3+SP3)*CP3 55 

CP5=(CP3+SP3)*(CP3-SP3) 56 

P5 1 = P2 1 *P4 1 — 0 .25714285*P31 * 57. 

DP 5 1 -P2 1 *DP41+QP2 1*P41-0.25714285*DP31 58 

P52=P21*P42-C.22857142*P32 59 

DP52 = P2 1 *DP42 + DP2 1 *P42-0 . 228571 42 »DP 32 60 

P53=P21*P43-0.14285714*P33 61 
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0P5 3 = P21 *DP43+DP2 1*P43— 0 . 1 42 85 714*DP33 62 

P54=P21*P44 63 

DP54=P21*DP44+DP21 *P44 64 

P55= r P22*P44 65 

DP55=4.G*P54 66 

AOR=AQR*AR 67 

C2=G(5»2)*CP2+H(5»2)*SP2 68 

C3=G( 5, 3) *CP3+H( 5,3 )*SP3 69 

C4=G(5,4)*CP4+H(5,4)*SP4 70 

C5 = G( 5» 5 ) *CP5 + H ( 5*5) *SP5 71 

BR=BR-5.0*A0R*(G( 5, 1 ) *P5 1 +C2*P 52+C 3*P5 3+C4*P 54+C5 *P55 ) 72 

BT=BT+A0R*(6< 5, 1 ) *0P5 1 +C2*DP 52 +C 3*0P5 3+C4*DP54+C 5*DP55 ) 7 3 

BP=BP-AGR*((G( 5,2)* SP2-HI 5,2)* CP2)*P 52+2.0* (G( 5,3)* SP3-HI 5,3) *CP3) 74 

1*P53+3.0*(G(5,4)*SP4-H(5,4)*CP4)*P54+4.0*<G(5,5)*SP5-H(5,5)*CP5)*P 75 

255) 76 

IF(NMAX- 5) 18,18,5 77 

C N= 6 78 

5 SP6=SP2*CP5+CP2*SP5 79 

CP6=CP2*CP5-SP2*SP5 80 

P61=P21*P51-0.25396825*P4l 81 

DP61=P21*DP51+DP21*P51-0.25396825*0P41 82 

P62=P21*P52-0.23809523*P42 83 

DP62=P2l*DP52+DP21*P52-0.23809523*DP42 84 

P63=P2l*P53-0. 19047619*P43 85 

DP63=P21*DP53+DP21*P53— 0.1904761 9 *DP43 86 

P64 = P21* P54— 0 .11 111111*P44 87 

DP64=p21*DP54 + DP21«P54-0. Ill 1 11U*DP44 88 

P65=P21*P55 89 

DP65=P21*DP55+DP21*P55 90 

P66 =P22* P 55 91 

0P66=5.0*P65 92 

AUR= ADR * AR 93 

C2=G(6,2)*CP2+H(6,2)*SP2 94 

C3 = G( 6, 3 )*CP3+H( 6,3)*SP3 95 

C4=G( 6,4 ) *CP4 + H t 6,4)*SP4 96 

C5=6I6,5)*CP5+H(6,5)*SP5 97 

C6=G(6,6)*CP6+H(6,6)*SP6 9b 

BR=BR-6.0*A0R*l G( 6, 1 ) *P6 1+C ?*P62 +C 3*P 63+C 4* P64+C 5 *P65 +C6* P66 ) 99 

BT=Br+AOR*{G(6, 1 ) *DP6 1 +C2* DP62+C 3 *DP63+C4*DP64+C 5 *DP6 5+C6 *0P66 ) 100 

BP=BP-AUR*( (G(6,2)*SP2-H(6,2)*CP2)*P62+2.0*{G(6,3)*SP3-H<6,3)*CP3) 101 

1*P63+3.0*(G16,4)*SP4-H{6,4)*CP4)*P64+4.0*(G(6,5)*SP5-H(6,5)*CP5) *P 102 

265+5. 0*(G{6,6)*SP6-H(6,6)*CP6)*P66) 103 

IF ( NM AX— 6 ) 18,18,6 104 

C N= 7 105 

6 SP7= ( SP4+SP4 )*CP4 106 

CP7=(CP4+SP4)*{ C P4— SP4 ) 107 

P7l=P21*P61— 0.25252525*P51 103 

DP71=P21*0P61+DP21*P6 1-0.2525 2525* DP 51 109 

P7?=P21*P62-0.24242424*P52 110 

0P72=P21*DP62+DP21*P62-0. 24242424* DP 52 111 

P73-P21*P63— 0.21212121*P53 112 

DP73=P21*DP63+DP21*P63-0.21212121*DP53 113 

P74=P21»P64-0.16161616«P54 114 

DP74-P21*DP64+DP2i*P64-0.16161616*DP54 115 

P75=P21*P6 5- 0. 090909 09»P55 116 

DP73=P21*DP65+DP21*P65-0.09090909*OP55 117 

P76=P21*P66 118 

DP76=P21*0P66+DP2 1*P66 119 

P77 = P22*P66 120 

DP77=6.0*P7fa 121 

AOR = A(JR*AR 122 



Table IV (Cont) 


C 

7 


c 

e 


C2=G(7,2)«CP2+H{7,2)*SP2 
C3=G( 7,3)*CP3+H(7,3)*SP3 
C4=G( 7,4)*CP4+H(7,4)*SP4 
C5=G« 7,5)*CP5+H(7,5)*SP5 
C6=G( 7,6)*CP6+H(7,6)*SP6 
C7=GI 7,7)*CP7+H(7,7)*SP7 

BR=BR-7.0*AOR*(G(7, 1 ) *P71+C2*P72+C3*P73+C4*P74+C5*P75+C6*P76+C7*P7 
17) 

BT=BT+AUR*(G(7,1 )*DP7l+C2*DP72+C3*DP73+C4*DP74+C5*DP75+C6*DP76+C7* 
1DP77) 

BP=8P-AGR*< (G( 7,2)*SP2-H(7,2)*CP2)*P72+2.0*<G(7,3)*SP3-H<7,3)*CP3) 
1*P73+3.0*(G( 7,4)*SP4-H( 7,4 )*CP4)*P74+4.0*(G( 7,5)*SP5-H( 7,5)*CP5) *P 
275+5. 0*(G1 7, 6)*SP6— HI 7,6 ) *CP6 ) *P 76+6. 0* ( G ( 7 , 7 ) *SP 7-H ( 7,7)*CP7)*P77 
3) 

IF ( NMAX- 7) 18,18,7 

N= 8 

SP8=SP2*CP7+CP2*SP7 

CP8=CP2*CP7-SP2*SP7 

P81=P21*P7l-0.25174825*P6l 

DP81=P21*DP71+DP21*P71— 0.25174825+DP61 

P82=P21*P72-0.24475524*P62 

0P82=P21*DP72+0P2 l*P72— 0.24475524*0P62 

P83=P2i*P73-0.22377622*P63 

DP83=P2 1*DP73+ DP 21*P7 3— 0.22377 622*DP63 

P84=P21*P74-0. 18881118*P64 

DP84=P21 +DP74+DP2 1*P7 4-0. 18881 1 18*DP64 

P8 5=P2 1*P75-0.1 39860 13*P65 

DP85 =P2i*0P75+DP21*P 75-0. 1398601 3* DP65 

P86=P2 1*P76— 0 .07692 30 7*P 66 

DP86=P21*DP76+DP21*P76-0.C7692307*DP66 

P87=P21*P77 

DP8f=P2l*DP77+DP21*P77 
P88=P22*P77 
0P88=7.0*P87 
A0R= A0R*AR 

C2=G( 8,2 ) «CP2+H( 8, 2)*SP2 
C3=G{8,3)*CP3+H(8, 3)*SP3 
C4=G( 8,4 )*CP4+H(8,4)*SP4 
C5=G( 3,5 )*CP5 + H( 8,5 >*SP5 
C6=G( 8,6) *CP6+H( 8,6)*SP6 
C7=G(8, H*CP7 + H(8,7)*SP7 
C8=G(8,8)*CP3+H(8,8)*SP8 

9R=BR-8.0*AQR*(G( 8, 1 ) *P8 1+C 2 *P 82+C3*P83+C4*P 84+C5* P85 + C 6*P86 + C7 * P8 
1 7+C8 • P88 ) 

BT-=BT+ACR*( G( 8 , 1 ) *DP8 1 +C2*DP32+C 3*DP 83 + C4*DP 84+C 5*DP35 + C6*DP86+C7 * 
IDP87+C8+DP88 ) 

9P=BP-AOR*t (G(8,2)*SP2-H(8,2)*CP2)*P8?+2.0+{G(8,3)*SP3-H{ 8,3)*CP3) 
1*P83+3.0*(G( 8,4) • SP4— H ( 8,4)*CP4)*P84+4.0*(G( 8, 5) *SP5-H 18,5)*CP5)*P 
285 + 5. 0*(G( 8,6 ) *SP6-H( 8,6)*CP6)*P86 + 6.0*(G(8,7)*SP 7-H( 8,7)*CP7)*P87 
3+7.0*(G( 8,8) *SP8-H( 8, 8 )»CP8) *P88) 

IF ( NMAX- 8) 18,18,8 

N= 9 

SP9=C SP5+SP5)*CP5 
CP9=(CP5+SP5)*(CP5-SP5) 

P9 1 = P2 1 *P8 1-0. 251 2820 5*P 71 

DP9l=P21*DP81+DP21*P81-0.25128205*DP71 

P92=P21*P82-0.24615384*P72 

DP 92=P21*DP82+DP21*P82-0. 24615384* DP72 

P93 -P2 1«P83— 0.230 76923*P73 

DP93=P2l*DP83+DP21*P83-0.23076923*DP73 

P94=P2 1 *P84— 0 .205 123 20* P 74 


123 
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1 75 

176 

177 

178 
1 79 
18C 
181 
182 
183 
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C 

9 


DP96=P2 1 *DP86+DP2l«P86-0.20512820*0P76 186 

P95=P21*P85-0. 169230 76»P75 185 

0P95=P21*0P35+DP21*P85-0. 16923076*OP75 186 

P96=P2l*P86-0. 1?307692*P76 187 

0P96^P2 1*DP86+DP 21 *P 86-0. 12 307692 *0P 76 188 

P97=P21*PB7-0.06666666*P77 189 

l)P9 7 = P21*DP3 7+DP2l*P87-0.06666666«DP77 190 

P98=P21 *P88 191 

DP93=P21*DP38+DP21*P88 192 

P99= P22 * P88 193 

DP99 = 8 . 0 *P98 196 

A0R=\0R*AR 195 

C2=G( 9,2 ) *CP2 + H ( 9,2 )*SP2 196 

C3-G ( 9 » 3 ) *CP 3+H ( 9 , 3 ) * SP 3 197 

C6=G(9,6)*CP6+H(9,6)»SP6 198 

C5 = G( 9, 5 )*CP5 + H(9,6)*SP5 199 

C6=G(9,6)*CP6+H(9,6)*SP6 200 

C7=G(9, 7)»CP7+H(9,7)*SP7 201 

C8=G(9,8)»CP8+H(9,8)*SP8 202 

C 9 = G ( 9 » 9 )»CP9+H(9,9)*SP9 203 

BR = 3R-9.0«Aur<*( G( 9, L ) * P9 1 +C2 *P 92+C 3* P9 3+ C 6*P 96+C5* P95 +C 6*P96+C 7* P9 206 

17+C8*P98+C9*P99) 205 

Br = C r+AUR*(G( 9, 1 ) *DP91+C2 *Dp92+C3*DP93+C6*DP96+C5*DP95+C6*DP96+C7* 206 

10P97+C8*0P93+C9*DP99) 207 

BP= B P-AUR* ( (0(9, 2 )*SP2-H(9,2 ) * CP2 ) * P9 2 + 2 . 0* ( G ( 9 , 3 > * SP 3-H ( 9,3 >*CP3) 2 08 

1*P93+3.G*(G( 9,6)*SP6-H(9,6 ) * CP6 ) * P96+6 . 0* ( G ( 9 , 5 ) *SP5-H(9,5 )*CP5) *P 209 

295+5.0* ( G( 9,6)*SP6~H( 9 , 6 ) * C P6 ) *P96+6 . 0* ( G ( 9 , 7 ) * SP 7-H ( 9 , 7 ) *CP 7 ) *P97 210 

3+7.0* (G( 9, 8) *SP8-H(9,3)*CP3) *P 98+ 8 . 0* ( G ( 9 , 9 ) * SP9-H ( 9 , 9 ) *CP9 ) *P99 ) 211 

I F ( NMAX- 9) 18, 18,9 212 

N= 1 0 213 

SP10'SP2*CP9+CP2*SP9 216 

CP10=CP2*CP9-SP2*SP9 215 

P101=P21*P9l-G.25093039*P8l 216 

DP101=P21*0P91+DP21*P91-0.25098039*DP81 217 

P102-=P21*P92-0.26705882*P82 218 

DP102=P21*DP92 +0P 2l*P9 2-0. 26705882* OP 82 219 

P10 3 = P2 l*P93-0. 2 352961 l*P8 3 220 

0P101=P21*0P93+0P21*P93— 0.23529611*DP83 221 

P106=P2l*P96-0.21563627*P86 222 

DP106=P2L*DP96+DP21*P96-0.21568627*DP86 2 23 

P105=P21*P95-0.18823529*P85 226 

DP105=P21*DP95+DP2i*P95-0. 18823529*DP85 225 

P106=P21*P96-0.15296117*P86 226 

0P106 = P21*DP96 + 0P21*P96-0.15296U7*DP86 22 7 

P107~P2l*P97-0.10980392*P87 228 

DP107=P21*DP97+DP21*P97-0. 10980392*DP87 229 

P108=P2l*P98-0.05882352*P88 230 

DP108=P21*DP98+0P21*P98-0.05882352*DP88 231 

P109=P2 1*P99 232 

DP109=P2i*OP99+OP21*P99 233 

P1010=P22*P99 236 

DP1010=9.0*P109 235 

AOR=AOR»AR 236 

C2=G( 10,2)*CP2+H( 10,2)*SP2 237 

C3 = G( 10,3)*CP3+H( 10,3)*SP3 238 

C6=G( 10,6)*CP6+H( 10,6 )*SP6 239 

C5=G( 10,5)*CP5+H( 10,5)*SP5 260 

C6=G( 10,6)*CP6+H( 10,6)*SP6 261 

C7=G( 10,7 )*CP7+H( 10,7)*SP7 262 

C8-=G( 10,8)*CP8+H( 10,8)»SP8 263 

C9=G( 10,9) *CP9+H( 10,9)*SP9 266 



Table IV ( Cont) 


C10=G( 10, 10 )*CP10+H( 10, 10)*SP10 24 5 

3R=8R-10.0*A0R*1G( 10, 1 ) *P 10 1+C2*P 1 02+C3* PI 03+C4*P 1 04+C5* PI 05+C6*P 1 246 

106+C7*P107+C8*P108+C9*P109+C10*P1010) 247 

RT =BF +AOR* (G( 10, 1 ) *DP 1 01 +C2*DP 102+C 3*DP 1 03+C4*DP 104+C 5*DP105 +C6*DP 248 

1106+C7*DP107+C8*DP108+C9*DP109+C10*DP1010) 249 

BP=BP— AUR* ( (GC 10,2)*SP2-HC 1 0 , 2 ) *CP2 ) » P102+2 . 0* < G ( 10 , 3 ) * SP3-H { 1 0 , 3 ) 250 

1*CP3)*P103+3.0*{G< 10, 4)*SP4-H( 10, 4 ) *CP4 ) *P 104+4. 0* ( G < 1 0 , 5 ) *SP5-H ( l 251 

20.5) »CP5)*P105+5.0*(G( 10,6 )*SP6-H( 10,6) *CP6)*P106+6.0*(G ( 10,7)* SP7 252 

i-Hl 10,7) *CP7)*P1 07+7.0*1 G( 10 , 8 ) *SP8-H ( 10 , 8 ) *CP8 ) *P 1 08+8. 0* ( G ( 10, 9 ) 253 

4*SP9-H( 10,9) *CP9)*P 109+9. 0*1 G( 10, 10)*SP10-H( 10, 10 ) *CP10 ) *P1010 ) 254 

IF(NMAX-IO) 18,18,10 255 

C N=1 1 256 

10 SPil=(SP6+SP6)*CP6 257 

CPll=tCP6+SP6)*(CP6-SP6) 258 

Plll=P2i*P101-0.25077399*P91 259 

DP111=P2 1*DP10 1+DP2 1*P 10 1-0. 25077399*DP91 260 

PI 12=P2 1*P102— 0. 2476 7801 *P92 261 

DP112=P21*DP102+DP21*P102-0.24767801*DP92 262 

P113=P21*P103— 0.2 3839009*p9 3 263 

DP113 ; =P2 1*DP103+0P2 1*P 1 03—0 .23 839009* 0P93 264 

P l 14=P2 1 *P 10 4-0. 2229 1021 *P 94 265 

DPI 14=P2 1 *DP 1 04+DP21*P 104-0. 22291 02 1*DP94 266 

PI 1 5= P2 1*P 105— 0 .20 123839*P95 267 

DP115=P21*DP105+DP21*Pl05-0.20123839*DP95 268 

P 1 16 =P21*P 106—0. 17337461*P96 269 

DPI 16=P21*0P106+DP21*P 106-0. 17 337461*0P96 270 

P117=P21*P 107—0 . 1393188 8* P97 271 

DPI 1 7 =P21*DP107+DP21*P10 7-0. 13931888* DP 97 272 

P118=P21*P108— 0. 09907120*P 98 273 

DP 11 3=P21*DP108+DP21*P 108-0. 099071 20*DP98 274 

P 1 19=P2 1*P 109— 0 .05263157*P99 275 

DP119=P21*DP109+DP21*P109-0.052631 57*DP99 276 

P1110=P2i*P1010 277 

DP 1 1 1 0=P2 1 *DP 10 10+DP2 1*P 1010 278 

P1111=P22*P1010 279 

DP1111=10.0*P1110 280 

AOR= AOR* AR 281 

C2=G< 11,2)*CP2+H( 11, 2)*SP2 282 

C3=G( ll,3)*CP3+Hl 11,3)*SP3 283 

C4=G( 11,4)*CP4+H( 11,4)*SP4 284 

C5=G( 11,5)*CP5+H( 11,5 )*SP5 285 

C6=G( 11 ,6) *CP6+H( 11,6 )*SP6 286 

C7=G( 11,7)*CP7+H( il,7)*SP7 287 

C8=G(11,8)*CP8+H(11,8)*SP8 288 

C9=G ( 11,9) *CP9+H( li,9)*SP9 28° 

C10=G( 11, 10)*CP10+H( 11,10)*SP10 290 

Cll = G(il,ll)*CPU+H(ll,li)*SPll 291 

BK=BR-11.0*A0R*)Gl 11, l ) *P1 1 1+C 2*P1 12 +C3*P11 3+C4*P 1 14+C5*P 1 1 5+C6*Pl 292 

116+C7*P117+C8»P118+C9*P119+C10*P1 1 10+C 1 1 *P 1 1 1 1 ) 293 

BT=BT+AOR*(G( 11, 1 ) *DP1 1 1 +C2 *DP 1 12+C3*0P 1 1 3+C4*DPl 1 4+C 5*DP 1 15+C6*DP 294 

lil6+C7*DPil7+C8*DP118+C9*DPU9+C10*DP1110+Cil*DPllll ) 295 

8P = BP— ACR*( I G I 11,2 ) *S P2— H ( 11,2 )*CP2)*P112+2.0*(G( 11,3)* SP3-H (11,3) 296 

i*CP3 )*P1 13+3.0 *(G( 11,4)* SP4-H( 1 1 , 4 ) *CP4 ) *P 1 14+4 . 0* ( G ( 1 1 , 5 ) *SP5-H ( 1 297 

21.5) *CP5)*PU5 + 5.0* (G( 11,6) *SP6-H( 11,6) *CP6)*P1 16 + 6. 0*(G( 11,7)* SP 7 298 

3-H( 11,7) *CP7)*P11 7+7.0* (G( 1 1 , 8 ) * SP8-H ( 1 1 , 8 ) *CP8 ) *P1 1 8+8 . 0* l G ( 1 1 , 9 ) 299 

4*SP9-H( 11,9) *CP9)*P 119+9.0* (G( 11, 10) *SP10-H( 11, 1 0 ) *CP 10 ) *P1 1 10+ 10. 300 

50* (G( 11, 11 )*SP11-H( 11, 11)*CP11)*P1111 ) 301 

IF(NMAX-ll) 18,18,11 302 

C N=1 2 303 

11 SP12=SP2*CP11+CP2*SP11 304 

CP12=CP2*CP11-SP2*SP11 305 



Table IV (Cont) 


P121=P21*P11 1—0 .2 50 6265 6* PI 01 306 

DP121=P21*DP111+DP21*P11 1-0 . 2 5062 656* DP 101 307 

P122=P2L*P112-0.24812030*P102 308 

DP122=P21*DP112+0P21*P112-0.24812030*DP102 309 

P123 : =P2I*P11 3-0. 240 60 1 5 0*P 103 310 

DP123=P2 1*DP l 13+DP21»P 11 3-0.24060 150*DP103 311 

P 1 24=P2 1 *P 1 l 4-0. 2280701 7*P 104 312 

DP 124=P2 1 *DP 1 14+DP2 1 *P 114-0. 2280 701 7*DP 104 313 

P125=P2L*P115-0.21052631*P105 314 

DP125=P21*DP115+DP21*Pl 15-0.2 10526 31*DP 105 315 

P126=P21*Pl 16-0. 18796992*P106 316 

DP126=P21*DP116+DP21*P116-0. 1 8 796 9 92* DP 106 317 

P 12 7=P21*P1 17-0. 16040 100*P 107 318 

DP127=P2 L *0P 1 17 +DP2 1 * PI 17-0. 16040 100*DP 107 319 

PI 28=P2 1*P 1 18-0. 12781 954*P 108 320 

DP128=P21*DP118+0P21*P113-0. 12781954»DP108 321 

P 1 29= P2 1 *P1 19-0. 09022 556*P 109 322 

DP129=P21*DP119+DP21*P119-0.09022556*DP109 323 

P121J=P21*P1110-0.04761904*P1010 324 

DP1210 = P21»DP1110+DP2l»PU10-0.04761904*DP1010 32 5 

P1211=P21*Pllll 326 

DP12I1=P21*DP1111+DP21*P1111 327 

P1212=P22»P1111 328 

DP1212=11.0*P1211 329 

AOR=AOR*AR 330 

C2=G( 12,2)*CP2+H( 12»2 )*SP2 331 

C3- G ( 12»3)*CP3+H( 12»3)*SP3 332 

C4=G ( 1 2 , 4 ) *C P4+H ( 12»4)*SP4 333 

C5=G< 12,5)*CP3+H(12,5)*SP5 334 

C6=G( 12,6)*CP6+H( 12,6)*SP6 335 

C7=G(12,7)*CP7+H( 12,7)*SP7 336 

C8=G( L2,8)*CP8+H< 12,8)*SP8 337 

C9=G< 12,9)*CP9+h( 12,9)*SP9 338 

C 10=G ( 12*10) *CP 1 0+H ( 12* 10)*SP10 339 

Cl i = G( 12,11 )*CP11+H( 12, 11 )*SPU 340 

C 1 2 = G ( 12,12) *CP 1 2+H ( 12,12 )*SP12 341 

8R=8R— 12.0*A0R*(G1 12, 1 ) *P 1 2 1 + C 2*P 1 22 +C3* P 12 3+C4 *P 1 24 + C5* PI 2 5+C6*P 1 342 

126+C7*P127+C8*P128+C9*P129+C10*P1210+Cll*P121 1+C12*P1212 ) 343 

BT = BT +AOR* ( G ( 12, 1 ) *DP12 1 +C2 *DP 122+C 3*DP 12 3+C4*0P 1 24+C5*0P 1 25 + C6*DP 344 

1126+C7*0P127+C8*DP128+C9»DP129+C10*DP1210+Cll*DP12il+C12»0P1212) 345 

BP = BP-AOR* ( C G ( 12,2)*SP2-H( 12 , 2 ) «CP2 ) * P 1 22+2 . 0* ( G ( 12 , 3 ) » SP3-H ( 12 , 3 ) 346 

1*CP3)*P 123+3. 0*(G( 12,4) *SP4-H( 1 2 , 4 ) *CP 4 ) *P 1 24+4. 0* { G ( 12 , 5)»SP5-H(1 347 

22,5) *CP 5) *P 125+5. 0*(G( 1 2, 6 ) +SP6-H ( 1 2 , 6 ) *CP6 ) * P126+6. 0* 1 G ( 12 , 7 ) * $P 7 348 

3— H ( 12,7 )*CP7)*P127+7. 0* ( G( 12,8)*SP8-H< 1 2 , 8 ) *CP 8 ) « P 128+8 . 0* ( G < l 2 , 9 ) 349 

4 + SP9-H { 12,9)*CP9)*P129+9.0*(G( 12, 10)*SP10-H( 12, 1 0 ) «CP 10 ) *P1 2 10+ 10. 3 50 

50* (G1 12,li ) *S P 1 1 — H ( 12, 11)*CP ll)*P 1211 + 11. 0*(G( 1 2 , 1 2 ) * SP 12-H ( 12 , 1 2 ) 351 

6*C P 1 2 ) *P 12 12 ) 352 

I F ( N M AX-12 ) 18,18,12 353 

C N=l 3 354 

12 SP13=( SP7+SP7 )»CP7 355 

CPl3=(CP7+SP7)*lCP7— SP7) 356 

P131=P21*P121-0.25051759*P111 357 

DP131=P2l*DPl2l+DP2l*P121-0.2505l759*DPlll 358 

P l 32 = P21*P 122-0. 24844 72 0*P L 12 359 

DP132=P21* DP 122+DP2l*P122-0. 2484472 0*DP112 360 

P 1 33= P2 1*P 12 3-0. 24223602*P 113 361 

DP133=P21*DP123+DP21*P123-0.24223602*DP113 362 

P134=P21*P 124—0. 23188405*P114 363 

DP134=P21*DP124+DP21*P124-0.23188405»DP114 364 

P135=P21*P125-0.21739130*P115 365 

DPI 35=P2 1 *DP125+DP21*P 125-0. 21739 130*DP 115 366 



Table IV (Cont) 


P136=P21*P 126-0. 19875776*P 116 367 

DP136=P2l*DP126+DP21*P126-0 .19875 776* DPI 16 368 

P137=P21*P127-0. 17598 343 *P 117 369 

DP137=P2i* DP 127 *0P2i*P127-0. 17598343* DP 117 370 

PI 38=P2 1 *P 128— 0 .14906832*P118 371 

DP133=P2l*DP128+DP21*P128-0. 14906832»0P118 372 

P139=P21*P129-0.11801242*P119 373 

DP139=P21*DP129+DP21*P129-0.11801242*DP119 374 

P1310=P21*P1210-0.08281573*P1110 375 

DP1310=P21*DP1210+OP21*P1210-0.08281573*DP1110 376 

P1311=P21*P1211-0.04347826*P1111 377 

DP131l=P21*OP1211+DP2i«P1211-0.04347826*DPllll 378 

P1312=P21*P1212 379 

DP1312=P21*0P1212+DP21*P1212 380 

P1313=P22*P1212 381 

PP1313=12.0*P1312 382 

A0R=AQR* AR 383 

C2=G(13,2)*CP2+H(13,2)*SP2 384 

C3=G( 13,3)*CP3+H( 13, 3)*SP3 385 

C4=G( 13,4) *CP4+H( 13,4)*SP4 386 

C5=G( 13,5)*CP5+H( 13,5)*SP5 387 

C6=G( 13,6)*CP6+H( 13,6)*SP6 388 

C7=Gt 13,7)*CP7+H( 13,7)*SP7 389 

C8=G (13,8) *CP8+H( 13,8)*SP8 390 

C9=G( 13,9)*CP9+H( 13,9)*SP9 391 

C10=G( 13,10)«CP 10+H{ 13, 101* SP 10 392 

C11 = G(13,11)*CPU+H( 13,U)*SP11 393 

C12=G( 13,12 )*CP12+H( 13, 12)*SP12 394 

Cl 3=G( 13,13)*CP13+H< 13, 13)*SP13 395 

BR=BR-13.G*A0R*(G( 1 3 , 1 ) *P 1 31 +C2*P 1 32+C3*P 1 33+C4*P134+C5*P135+C6*P 1 396 

136+C7*P137+C8*P138+C9*P139+C1C*P1310+C11*PI 31 i+C 12*P1 312+C1 3*P1 31 3 39 7 

2) 399 

BT=BT+AUR*(G{ l 3, l ) *OP 1 3 1+C2*DP 1 32+C3*DP 133+C4*DP134+C5*DP135+C6*DP 399 

1136+C7*DP137+C8*DP138+C9*DP139+C10*DP1310+C11 *DPl 31 i+C12*DP 1 3 12+C 1 400 

23*DP 1 3 1 3 ) 401 

BP=BP-AOR*( (G( 13,2)*SP2-H< 1 3, 2 ) *CP2 ) *P1 32+2.0* (G( 13,3 )*SP3-H( 13,3 ) 402 

1*CP3 )*P1 33+3.0* (G( 13,4) * SP4-H ( 13, 4 ) *CP4> *P1 34+4 . 0* CG ( 1 3 , 5 ) *SP5-H ( 1 403 

23,5 )*CP5) *P1 35+5.0* (G( 13,6) *SP6-H( 1 3, 6 > *CP6 )*P1 36 + 6. 0* (G ( 13 , 7 I * SP7 404 

3— H (13,7)*CP7)*P137+7.0*(G( 1 3, 8 ) •SP8-H ( 1 3 , 8 ) *CP8 ) *P1 38+8.0* ( G J 1 3,9 ) 405 

4* SP9-H ( 13,9)*CP9) *P139+9.0*(G( 1 3, 10 )*SP 10-H1 1 3, 10 ) *CP10) *P1310+ 10. 406 

50* (G( 13, 11 )*SPil-H( 13, 11 )*CP11 ) *P 1311+11. 0*(G( 13, 12 ) * SP12-H ( 13, 12) 407 

6*CP1 2 ) *P 1312+12. 0* ( G ( 13, 13)*SP13-H< 13, 13 )*CP13)*P1313) 408 

If ( NM AX— 13 ) 18,18,13 409 

N=14 410 

SP14=SP2*CP13+CP2*SP13 411 

CP14=CP2*CP13-SP2*SP13 412 

P141=P21*P1 31-0. 25043478*P12l 413 

DP141=P21*DP131+DP21*P131-0.25043478*DP121 414 

P142=P21*P 132-0. 24869 565*P 122 415 

DP142^P21*DP132+DP21*P132-0.24869565*DP122 416 

P143-P2l*P 133—0 .24347826*P123 417 

0P143=P21*DP133+DP21*P133— 0.24347826*DP123 418 

P144=P21*P134-0.23478260*P124 419 

DP144= P21* DP 134+0P21*P134- 0.23478260* DP 124 420 

P145=P21*P135-0.22260869*P125 421 

DP145=P21*DP135+0P21*P135-0.22260869*DP125 422 

P146=P21*P136-0.20695652*P126 423 

DP146=P2l*DP136+DP21*P136-0. 20695652*DP126 424 

P147=P21*P137-0. 18782 60 8*P127 425 

DP147=P21*DP137+DP21*P137— 0.18782 60 8*DP127 426 

PI 43^ P2 1 *P 1 38-0. 16521 739*P 128 427 



Table IV (Cont) 


DP143=P21*nP138+DP2i*P138-0. 1 652 I 739* DP 1 2 8 42 8 

P149=P21*P139— 0.13913043*P129 429 

DP14^ = P21*DPl39+'DP21*P139-0. 1391 304 3*DP129 4 30 

P1410=P21*P1310-0.10956521*P1210 431 

nP1410=P21*DPl 310+UP21*P 1310-0. 10956521 *DP12 10 432 

Pl41l=P2l*P13ll-0. 076521 73*P1211 433 

DP 14 1 1=P21*DP1311+DP21*P1311-C.07652173»DP1211 434 

P1412=P21*P1312-0.04000000*P1212 435 

DP14l2=P21*0Pi312+DP2i*Pl312-0.04000000*DP1212 436 

P141 3=P21*P1313 437 

DP14l3=P21*DP1313+0P?l*P13l 3 438 

P1414=P22*P1313 439 

DP1414=13.0*P1413 440 

A0R= AOR * AR 441 

C2-G ( 14*2)*CP2+H( 14»2)*SP2 442 

C3=G( 14,3)*CP3+H( 14, 3)*SP3 443 

C4=G(14,4)*CP4+H( 14,4)*SP4 444 

C5 = G< 14,5 )*CP5+Hl 14,5 )*SP5 445 

C6=G( 14,6)*CP6+H( 14,6)*SP6 446 

C 7=G ( 1 4 , 7 ) *CP 7+H { 14 » 7 ) *$P 7 44 7 

C8 = G( 14,8)*CP9+H( 14*8)*SP8 448 

C 9 = G t 14,9)*CP9 + H( 14,9)*SP9 449 

C 10 = G( 14, 10) *CP10+H( 14,10) *SP 10 450 

C11 = G{ 14, 11 )*CP11+H( 14, 1 1 ) * SP 1 1 451 

C12=G( 14,12)»CP12+H( 14, 12)*SP12 452 

C 1 3=G ( 14,13) *CP 1 3 +H( 14,13)*SP13 453 

C14=G(14,14)*CP14+H( 14, 14)*SP14 +54 

9R=BR— 14.G*AOR*(G( 14, 1 ) *P 14 1+ C 2*P 1 42 +C 3 *P 143+C 4*P 144+C5 *P 14 5+C 6*P 1 455 

146+C 7 *P 147+C8*P 1 48+C9*P149+C 10*P 1 4 1 0+C 1 1 *P14 ll+C 1 2*P1 41 2 + C l 3 *P 14 1 3 456 

2+C 14»P1414) 457 

BT = R T + AGR* ( G 1 14, 1 ) *DP 14 1 +C2 *DP 14 2+C 3*DP 143 + C4*DP 1 44+C 5*DP 145 + C6*DP 45 8 

1 1 46+C 7»DP 147 <-C8*0P 1 48+C9*DPl 49+C 1 0*DP 14 10+C 1 1 *DP 1 41 1 + C 1 2*DP 1 412+C 1 459 

23*DP1413+C14*QP1414) 460 

EP=nP-AOR*( ( G ( 14, 2 )*SP2-H( 14,2 ) *C P 2 ) *P 142 +2 . 0* l G ( 1 4 , 3 ) * SP3-H ( 14,3) 461 

1*CP3 ) *P 143+3.0* (G( 14,4) *SP4-H( 1 4, 4 ) »C P4 ) * P 1 4 4+4. 0* ( G ( 14 , 5 ) * SP5-H ( 1 462 

24,5) *CP5)*P145+5.0*(G( 14,6)*SP6-H( 1 4 , 6 ) * C P6 ) *P 1 46+6 . 0* ( G ( 1 4 , 7 ) * S P 7 463 

3-H( 14 , 7 ) *CP 7 ) *P 14 7 + 7 . 0* ( G ( 1 4, 8 ) *SP8— H ( 1 4 , 8 ) *CP8 ) * P 148 + 8 . 0* ( G ( 14,9) 464 

4*SP9-H( 14,9)*CP9)*P149+9.0*(G( 14, 10)*SP10-H( 14 , I 0 ) *CP 10 ) * PI 4 1 0+ 10 . 465 

50* (G( 14, 11 )*SP11-H( 14, 11 )*CP11 ) *P 141 1 + 1 1.0* (G( 14, 12 )*SP12-H< 14,12 ) 466 

6*CP12) *P1412+12.0*(G( 14,13 ) *SP 1 3-H ( 14 , 1 3 ) *CP 13 ) *P 141 3+1 3. 0* ( G ( 14 , 1 467 

74 ) *5P 1 4-H( 14, 14)*CP14)*P1414) 466 

IFCNMAX-14) 18,18,14 469 

C N=l 5 470 

14 SP 15= ( SP8+SP8 ) *CP9 471 

CP15= (CP8+SP8 )*( CP8-SP8) 472 

P151=P21*P141— 0.25037037*P131 473 

DP151=P21*DP14i+0P21*P141-0.25037037*DP131 474 

P152=P21*P 142-0. 2488 8888 *P 132 475 

DP 1 52 = P2 1 *DP 142 + DP2 1 *P 142-0. 24 888888*DP 132 476 

P 1 5 3 = P2 1*P 14 3-0. 2 4444444* PI 33 47 7 

DP153=P2i*DP143+DP21»P143-0.24444444*0Pl 33 478 

P154=P21*P144-0.23703703*P 134 479 

DP154=P21*DP144+0P2l*P 144-0. 2370 3703*DP 1 34 480 

P155 = 1>21*P 145-0. 2 2666666 *P 135 481 

DP155 = P21* DP 145 + DP21*P14 5- 0.22666666* DP 135 482 

P 1 56 = P2 l *P 146-0. 2 1333333* PI 36 48 3 

DP156-P21*DPl46+DP21*P146-0.21333333*DP136 484 

P 157=P21*P 147-0. 19703703* P 137 485 

OP157=P21*DP147+DP21*P147-0. I 97037 03*DP 1 37 486 

PI 58=P2 1 *P 148-0. 177 777 77*P 138 487 

DP153=P21*DP148+DP21*P148-0. 1 7777777*DP 1 38 488 


I 


Table IV (Cont) 


P159=P2i*P149— 0. 15555555*P139 489 

DP159=P2 1*DP 149+DP21*P 149— 0. I5555555*0P 139 490 

P1510=P2 1*P1410— 0 . 130 3/03 7»P1310 491 

DP1510=P2 i*DP 1410+DP2 1*P1410— 0 . 1 3037037*DP1310 49 2 

P15 11=P21*P141 1-0. 10222222 *P 1311 493 

OP 151 l=P21*DP1411+DP21»P14il— 0.10222222*DP1311 494 

PI 512=P2 1*P1412-0.071 lllll*P1312 495 

DP1512=P21*DP1412+DP21*Pl412-0.07illlll*DP1312 496 

P1513=P2i#P1413-0.03703703*P1313 497 

DP 151 3=P21*DP1413+0P2 1*P 1413-0. 03703703*DP1313 498 

P1514=P21*P1414 499 

DP1514=P21*DP1414+DP21*P1414 500 

P1515=P22*P1414 501 

DP1515=14.0*P1514 502 

A0R= A0R*AR 503 

C2=G( 15,2)*CP2+H( 15,2)*SP? 504 

C3=G( 15,3)*CP3+H(15,3)»SP3 505 

C4=G( 15,4)*CP4+H( 15,4)*SP4 506 

C5=G( 15,5)*CP5+H( 15,5)*SP5 507 

C6=G( 15,6)*CP6+H( 15,6)*SP6 508 

C7=G( 15,7)*CP7+H< 15,7 )*SP7 509 

C8=G( 15,3)*CP8+H( 15,8)*SP8 510 

C9=G( 15,9)*CP9+H< 15,9)*SP9 511 

C10=G( 15, 10)*CP10+H< 15, 10)* SP 10 512 

C11=G< 15,11)*CP11+H115,U)*SP11 513 

C12=G( 15,12 )*CP12+H< 15,12)*SP12 514 

C13=G( 15, 13 ) *CP 1 3+H( 15, 13)*SP13 515 

C14=Gl 15,14) *CP14+H( 15, 14)*SP14 516 

C15=G( 15,15) *CP 15+H( 15, 15)»SP15 517 

BR=BR— 15.0*A0R*( G( 15 , 1 ) *P 1 5 1+C2*P 1 52+C3*P 153+C4*Pi 54+C5 *P 155+C6* PI 518 

156+C7*P157+C8*P158+C9*P159+C10*P1510+C11*P1511+C12*P1512+C13*P1513 519 

2+Cl4*P1514+C15*P1515) 520 

Br=8T+AOR*(G( 15, 1 )*DP15l+C2*DP152+C3*0P153+C4*DP154+C5*DP155+C6*DP 521 

1156+C7*0P157+C8*DP158+C9*DP159+C10*DP1510+C1 1*DP1 51 1+C12*DP1512+C1 522 

23*DP1513+C14*DP1514+C15*DP1515) 523 

BP=BP— AOR* ( (G( 15,2 )*SP2-H( 15,2 ) *CP2 > * PI 52+2 . 0* ( G C 15 , 3 ) *SP 3-H { 15 , 3 ) 524 

1*CP3)*P1 53+3.0 *<G( 15,4 )*SP4-H( 15, 4 ) *CP4 ) *P1 54+4. C* ( G ( 1 5 , 5 ) *SP5-H ( 1 525 

25,5) *CP5)*P155+5.0*(G< 15,6) *SP6-H( 15,6 )*CP6)*P156+6.0*(G( 15,7) *SP 7 526 

3-H (15,7)*CP7)*P157+7.0*(G( 1 5 , 8 ) *SP 8-H ( 15 , 8 ) *CP8 ) *P 158+8 . 0* ( G ( 15 , 9 ) 527 

4»SP9— H( 15,9) *CP9 ) *P 15 9+9. 0* ( G( 15, 1Q)*SP10-H( 15, 10 )*CP10 ) *P15 10+ 10. 528 

50* (G( 15, 11 )*SP11-H( 15, 11 )*CPU)*P1511 + li.0*(G< 15 , 12 )*SP l?-H < 15,12) 529 

6*CP12)*P1512+12.0*(G( 15,13)*SP13-H( 15, 13 ) *CP l 3 ) *P 15 1 3+1 3. 0* ( G < 1 5 , 1 530 

74) *SP14-H( 15, 14) *CP 14 )*P1 514+14.0* (G( 15, 15) *SP15-H( 1 5 , 1 5 ) *CP 1 5 ) *P 1 531 

8515) 532 

I F ( MM AX— 15 ) 18,18,15 533 

C N=16 534 

15 SP 16=SP2*CP15+CP2*SP15 535 

CP16=CP2*CP15— SP2*SP15 536 

P161=P21*P151-0.25031928*P141 537 

DP161=P?1* DP 151+DP2i*P151-0. 2503192 8* DP141 538 

P 16 2= P21 *P1 52-0. 249042 14*P 142 539 

DP 1 6?=P2 1 *DP 152+DP21*P 152-0. 249042 14*DP 142 540 

P163=P21*P 153-0.2452 1072*P143 541 

DP16 3 = P21*DP1 53+DP2 1* P 15 3-0 .24521072* DP 143 542 

P164=P21*Pl 54-0 . 23882503*P144 543 

DP164=P21*DP154+DP21*P154-0.23882503*DP144 544 

P 1 65=P2 1»P 155-0. 22988505*P 145 545 

DP165=P21*DP155+DP21*P1 55-0 .22988505* DP 145 546 

P166=P21*P 156—0. 21839080*P 146 547 

DP166=P21*DP156+DP2 1*P 156-0. 21 839080* OP 146 548 

P167 = P21*P15 7- 0.20434 2 27*P 147 549 


Table IV (Cont) 


DP167=P21*DP157+DP21*P157-0.20434227*DP147 550 

P 1 68=P2 1*P 158-0. 18773946*P 148 551 

DP168=P21*DP 1 58+DP2 1*P158-0.18773946*DP148 552 

P 169= P2 1 *P 159-0. 16858237*P 149 553 

DPI 69=P21*DP159+DP21*P1 59-0. 16858237*DP 149 554 

PI 610=P2 1*P1 5 10-0. 146871 00*P 1410 555 

DP1610=P2l*DP 15 10+DP21*P 1510-0. 14687 100«DP 1410 556 

P1611=P21*P1511— 0.122 60 536*P1411 557 

DP161 i=P21*DP1511+DP21*P 15 1 1-0. 12260536*DP141 1 558 

P1612=P21*P1512-0.09578544*P1412 559 

DP1612=P2i*DP1512+DP21*P1512-0 . 09578544*DP14 12 560 

P1613=P2l*P1513-0.06641123*P14l3 561 

DPI 61 3=P2 1*DP 15 13+DP2i*P 15 13-0.066411 2 3*DP14l 3 562 

PI 614=P21*P1 5 14-0.03448 275*P1414 563 

DP 16 14= P2 1 *DP 15 14 +DP2 1*P 15 14-0. 03448 2 75*DP 14 14 564 

P1615=P21*P1515 565 

DP1615=P21*DP1515+DP21*P1515 566 

P1616=P22*P1515 567 

DP1616=15.0*P1615 568 

A0R= AOR* AR 569 

C2=G( 16,2)*CP2+H< 16,2 )*SP2 570 

C3=G( 16,3)*CP3+H{ 16,3)*SP3 571 

C4=G(16,4)*CP4+Hl 16,4)*SP4 572 

C5=G( 16,5)*CP5+H< 16,5 )*SP5 573 

C6=G( 16,6)*CP6+H( 16,6)*SP6 574 

C7=G( 16, 7)*CP7+H( 16,7)*SP7 575 

C8=G( 16,8)*CP8+H( 16, 8 )*SP8 576 

C9=G« 16,9)*CP9+H( 16,9)*SP9 577 

C10=G( 16,10 ) *CP 10+H ( 16,10) *SP10 578 

Cll=G< 16,11 )*CP11+H( 16,11)*SP11 579 

C12=G(16,12)*CP12+H( 16, 12)*SP12 580 

C l 3=G ( 16,13) *CP 1 3+H( 16,13)*SP13 581 

C14=G(16,14)*CP14+Hl 16,14)*SP14 582 

C 1 5 = G ( 16, 15 )*CP15+H( 16, 15)*SP15 583 

C 1 6=G ( 16,16) *CP16+Hi 16, 16)* SP 16 584 

BR=BR-16.0*A0R*(G( 16, 1 ) *P1 61 +C2*P 162+C 3*P 163+C4* P 164+C 5*P 1 65+C6* P 1 585 

166+C7*P167+C8*P 168+C9*P 169'+C 10*P 1 61 0+C 11*P1611+C12*P1612+C13*P1613 5 86 

2+C 1 4*P 1 6 14+C 15*P1615+C16*P1616) 587 

BT=BT+ACR*(G(16, 1 ) *DP 16 1+C2 *DP 162 +C 3*DP 163+C4*DP164+C 5*DP 1 65 +C6*DP 588 

1166+C7*DP167+C8*DP168+C9*DP169+C10*DP 1610+C11*DP161 1+C12*DP 161 2+C l 589 

23*DP1613+C14*DP 16 14+C 15*DP1615+C16*DP1616) 590 

BP=BP-AOR*( ( G ( 16,2)*SP2-H( 16,2)*CP2 ) *P 162+2. 0* ( G ( 16 , 3 ) * SP 3-H ( 1 6 , 3 ) 591 

l*CP3)*P163+3.0*(G( 16,4)*SP4-H( 16, 4 ) *C P4 ) *P 164+4. 0* ( G ( 1 6 , 5 ) *SP5-H ( 1 592 

26, 5) *CP5)*P165 + 5.0*(G( 16,6 )*SP6-H{ 1 6 , 6 ) *CP6 ) * P 166 + 6. 0* ( G ( 16,7)*SP7 593 

3-H( 16,7 )*CP7)*P1 67+7.0* (G( 16 , 8 ) *SP 8-H ( 16 , 8 ) *CP8 ) *P 168+8 . 0* ( G l 16 , 9 ) 594 

4*SP9-H( 16,9) *CP9) *P 169+9. 0*1 G( 16, 10 )*SP10-H( 16, 10 ) *CP10 ) *P1 61 0+ 10 . 595 

50* (G( 16, 11 )*SP1 1— H ( 16, 11 )*CP11 )*P1611+11.0*(G( 16, 12)*SP12-H( 16,12) 596 

6*CP 12 )*P1612+12.0*(G( 16,13 )*SP13-H{ 16,13 )*CP13)*P1613+13.0*(G( 16, 1 597 

74) *SP14-H( 16, 14)*CP14)»P1614+14.0*(G( 16 , 1 5 ) * SP 1 5-H ( 16 , 1 5 ) *CP 1 5 ) *P1 598 

8615+15. 0*(G(16, 16 )*SP16-H( 16, 1 6 ) *CP 16 ) *P 1 61 6 ) 599 

IFINMAX-16) 18,18,16 600 

C N= 17 601 

16 SP17=(SP9+SP9)*CP9 602 

CP17=(CP9+SP9 )*ICP9-SP9) 603 

P171=P21*P 161-0. 25027808*P151 604 

DP171=P21*DP16l+DP2l*P161-0.25027808*DP151 605 

P172=P21*P162-0.24916573*P152 606 

DP172=P2l*DP162+DP21»P162-0.24916573*DP152 607 

P 1 73=P21*P 16 3-0. 24582 869*P 153 608 

DP173=P21*DP163+DP21*P163-0.24582869*DP153 609 

P174=P21*Pl64 -0.2402669 6* P154 610 


1*4 >u *■“* c\i ^ *«T 


Table VI (Cont) 


i C EPCCb= 19A6.0 


VESTINE 


N h G H G T 


HT GTT 


2 1-3C57C.CCC0 

2 2 -21 1C * CCCG 

3 1 -1270. COCO 

3 2 2960. CCCC 

3 2 164C.CCCC 

4 1 115C.CCCC 
A 2 -173C.CCC0 

4 3 121C.CCC0 

A 4 68C.CCCC 

5 1 920. COCO 

5 2 780.CCC0 

5 3 580.CCCC 

5 4 — 38C .COCO 

5 5 310. COCO 

6 1 -240. CCCC 

280. CCCC 
2CC.CCCC 
-6C.CCCC 

-15C.CCCG 
-7C.CCCC 
6C.CCC0 
110.CCCC 
-3C.CCC0 
-26C.CCC0 
-3C.CCC0 

7 6 30.CCC0 

7 7 -11C.CCCC 


-0. 

5.CCCC 

-0. 

-c. 

5810. OCCO 

2. COCO 

l.CCOO 

-c. 

-0. 

-18.CCCC 

-c. 

-c. 

166C.0CC0 

c. 

-20. CCCC 

— ■ c * 

5AC.CCC0 

2. COCO 

-14.CCCC 

-o. 

-c. 

-c. 

-0. 

-0. 

-52C.CCCO 

-c. 

-c. 

— c • 

19C.CCC0 

-c. 

-c. 

-c. 

3C.COCO 

-c. 

-c. 

-0. 

-0. 

-c. 

-0. 

-0. 

1A0.CCC0 

-c. 

-0. 

-0. 

-28C.C0CC 

-c. 

-c. 

-c. 

-80. OCCO 

-0. 

-c. 

-o. 

-130. COCO 

-c. 

-0. 

-0. 

-0. 

-0. 

-c. 

— c . 

30. COCO 

-c. 

-c. 

-0. 

7C.CCC0 

-c. 

-c. 

-0. 

-20. COCO 

-c. 

— c • 

-0. 

-1AC.CCCO 

-c. 

-c. 

-0. 

8C.CCCO 

-c. 

-c. 

-c. 

-0. 

-c. 

— c • 

-0. 

-70. COCO 

-0. 

-c. 

-0. 

150. COCO 

-0. 

-c. 

-c. 

-20. COCO 

-c. 

-c. 

-o. 

0. 

-G. 

-c. 

-0. 

- 1 C. coco 

-0. 

-c. 

-0. 

-AO. COCO 

-c. 

-c. 

-0. 


HTT 


Ve stine (I960) Table 1 correcting G(7,6) from -130. to +30 
from Table 3 for 1942.5. 


Secular change 


ooooooooooooooooooooooooooo 


Table VI (Cont) 


1 

C 

EPCCH= 

1958. 5 

NA GATA 

CGCT I 

( INTERNAL ) 


N 

h 

G 

H 

GT 

HT 

GTT 

HIT 

2 

1 - 

3C45C.CCCC 

- 0 . 

13.6CCC 

-C. 

-C. 

-c. 

2 

2 

-222C.CCC0 

5840. CCCC 

5.5CCC 

2.3CCC 

-C. 

-c. 

3 

1 

- 1 5 1C . CCCC 

-c. 

-21.5CCC 

— C . 

- 0 . 

-c. 

3 

2 

295C.CCC0 

— 194C.0CGC 

-i.eocc 

-16. CCCC 

- 0 . 

-c. 

3 

3 

149C.C0C0 

210. CCCC 

“C. 1CCC 

-17. 1CCC 

- 0 . 

-c. 

4 

1 

116C.CCCC 

- 0 . 

1.9CCC 

- 0 . 

-C. 

-c. 

4 

2 

— 166C . CCCC 

-45C.CCC0 

-R.2CCC 

6.2CCC 

- 0 . 

-c. 

4 

3 

1230. COCO 

240.CCC0 

1.1CCC 

3.5CCC 

- 0 . 

-c. 

4 

4 

930. CCCC 

-9C.0CCC 

-0.5CCC 

-6.7CC0 

- 0 . 

-c. 

5 

1 

94C.CCC0 

-C. 

-A.4CCC 

-c. 

-C. 

-c. 

5 

2 

750. CCCC 

190.CCC0 

2.7CCC 

C. 

-c. 

-c. 

5 

-a 

64C.CCC0 

“230 . OCCO 

-2. CCCC 

-1. 1 CCC 

-c. 

-c. 

5 

4 

-37C.CCC0 

-30. CCCC 

-l.CCCC 

3.5CCC 

- 0 . 

-c. 

5 

c. 

25C.CCCC 

-2C0.CCC0 

-2.2CCC 

-2. 1CCC 

-c. 

-c. 

6 

1 

-19C.CCCC 

- 0 . 

2.1CCC 

-C. 

-C. 

-c. 

6 

2 

34C.CCC0 

10. COCO 

—l.CCCC 

i.eccc 

-c. 

-c. 

6 

3 

210.CCC0 

90. OCCO 

C.2CCC 

1.7CC0 

- 0 . 

-c. 

6 

4 

-80.CCC0 

-9C.CCC0 

l.CCCC 

-3.8CCC 

-G. 

-c. 

6 

c 

- 1 5C • CCCC 

-120. CCCC 

-C.5CCC 

C.7CCC 

- 0 . 

-c. 

6 

6 

-90. CCCC 

ICC. CCCC 

-C.2CCC 

0.2CCC 

-c. 

- 0 . 

7 

1 

70.CCC0 

- 0 . 

1.7CCC 

-C. 

- 0 . 

-c. 

7 

2 

4C.CCC0 

-5C.CCC0 

2.1CCC 

-C.7CCC 

-G. 

-c. 

7 

3 

-10.CCCO 

60. CCCC 

C.1CCC 

-C.2CCC 

- 0 . 

-c. 

7 

4 

-200.CCC0 

20.CCC0 

-C.2CCC 

1.5CCC 

- 0 . 

- 0 . 

7 

c 

-50.CCC0 

-10.CCC0 

1.50CC 

-1. 1CC0 

-c. 

-c. 

7 

e 

-2C.CCC0 

-1C.CCC0 

C.5CCC 

O.7CC0 

-c. 

-c. 

7 

7 

-ICC. CCCC 

-I0.CCC0 

-C. 6 CCC 

l.CCCC 

-c. 

-c. 


Nagata and Oguti (1962) 



Table VI (Cont) 


I 

C 

E PCCH = 

1965. C 

LEATCN 

EVANS (X+Y> 



N 

H 

G 

h 

GT 

HI 

GT T 

HTT 

2 

1 - 

-30375. COCO 

- 0 . 

15.5CCC 

-C. 

-c. 

-c. 

2 

2 

-2087. CCCO 

5769. CCOO 

8.30CC 

0.6CG0 

— C • 

-c. 

3 

1 

— 1648.CCC0 

- 0 . 

-26.6CCC 

-C. 

-c. 

-c. 

3 

2 

2954. CCCC 

-1995. CCCO 

- 1.3CCC 

-1 1.4CCC 

-c. 

-c. 

3 

? 

1579. CCCO 

116. COCO 

1.3CCC 

-18.2CCC 

- 0 . 

-c. 

4 

1 

1164. CCCO 

- 0 . 

-c. 

-c. 

- 0 . 

-c. 

4 

2 

-2C33.CCCC 

-389. CCCO 

-9.5CCC 

3.2C0C 

-c. 

-c. 

4 

2 

1295. COCO 

230. CCCO 

-1.9CCC 

1.6CCC 

-c. 

-c. 

4 

A 

8 80. CCCC 

-141. COCO 

-C. 6 CCC 

-8.5C00 

-c. 

-c. 

5 

1 

93C.CCCG 

- 0 . 

C. 6 CCC 

-C. 

- 0 . 

-c. 

5 

2 

811. CCCC 

142. CCCO 

l.CCCC 

3 • C CCG 

-c. 

-c. 

5 

3 

45C.CCC0 

-276.CC0C 

-2 . 2CCC 

-C. 7CCC 

-c. 

— c • 

5 

A 

-402. CCCO 

5. COCO 

C.20CC 

2. 7CC0 

- 0 . 

-c. 

5 

c 

262. CCCO 

-264. CCCO 

-3. CCCC 

-2.7CCC 

- 0 . 

-c. 

6 

1 

— 175. CCCC 

-C. 

C . 8 CC C 

-c. 

~e. 

-c . 

6 

2 

357. CCCC 

30. CCCC 

C.4CCC 

1.9CCC 

-c. 

-c. 

6 

0 

248. CCCO 

135. CCCO 

1.6CCC 

2. 3CCG 

- 0 . 

-c. 

6 

A 

-2C.CCCC 

-123. CCCO 

-C.30CC 

-1.8CCC 

-c. 

-c. 

6 

c 

-171. CCCO 

-ICC. COCO 

-1.1CCC 

1.4CCC 

-c. 

-c. 

6 

6 

-64. CCCO 

84. CCCC 

1.7CCC 

C.5CCC 

- 0 . 

-c. 

7 

1 

42. CCCO 

- 0 . 

-C. 

-c. 

- 0 . 

-c. 

7 

2 

55. CCCC 

-18. COCO 

-C.5CCC 

-2.2C0C 

-c. 

-c. 

7 

3 

12. CCCC 

101. CCCO 

1.8CCC 

C.2C0C 

-c. 

-c. 

7 

A 

-239. CCCO 

60. CCCO 

1.4CCC 

C.9CCC 

- 0 . 

-c. 

7 

5 

16. CCCO 

-32.0000 

C. 6 CCC 

-2. 1CC0 

- 0 . 

-c. 

7 

6 

6 . CCCO 

-27. CCCO 

C.2CCC 

-0. 1 CCG 

-c. 

-c. 

7 

7 

-lie. CCCO 

-12. CCCO 

-2.2CCC 

-C.2CCC 

-c. 

-c. 

8 

i 

77. CCCO 

- 0 . 

-C. 

- 0 . 

- 0 . 

-c. 

8 

2 

-56. CCCO 

-47. CCCO 

-c. 

-c. 

- 0 . 

-c. 

8 

3 

8 . CCCO 

-35. CCCC 

- 0 . 

-c. 

-c. 

-c. 

8 

A 

5. CCCO 

-9. CCCC 

-c. 

-c. 

-c. 

-c. 

8 

5 

-35. CCCC 

2.0CC0 

-c. 

-c. 

- 0 . 

-c. 

8 

6 

-16. CCCO 

27. CCCC 

-c. 

-c. 

- 0 . 

-c. 

8 

7 

9. CCCC 

-17. COCO 

-c. 

-c. 

-c. 

-c. 

8 

e 

2. CCCO 

-24. COCO 

-C. 

-c. 

-c. 

-c. 

9 

1 

11. CCCO 

- 0 . 

-C. 

-G. 

- 0 . 

-c. 

9 

2 

23.0CCO 

5.C0CC 

-c. 

-c. 

-c. 

-c. 

9 

3 

- 6 . CCCC 

-7. CCCC 

- 0 . 

-c. 

-c. 

-c. 

9 

A 

-17. CCCO 

3. CCCO 

-C. 

-0 • 

— c . 

-c. 

9 

G 

5. CCCO 

-2C.00C0 

-C. 

-c. 

- 0 . 

-c. 

9 

6 

2C.CCC0 

8 . CCCO 

-c. 

-c. 

-c. 

-c. 

9 

7 

-l.CCCC 

26. CCCO 

-C. 

-c. 

-c. 

-c. 

9 

8 

12. CCCC 

1C. CCCC 

-C. 

-c. 

-c. 

-c. 

9 

9 

7. CCCC 

-12. CCCC 

-C. 

-c. 

-0. 

— c • 


Private Communication (1964) 
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Table VI (Cont) 


C 1 E PCCH= 1942.5 JCNES NELCTTE ( OBLATE ) 

N fc C H GT HT GTT HTT 

2 1 3022C.CCCC -C. -C. -C. -C. -C. 

2 2 21bC.CCCO -5 53C . CCCC -C. -C. -0. -C. 

3 1 I76C.CCCC -0. -C. -C. -C* -C . 

-5C7C.CCCC 2590. OCCC' -C. -C. -0. -C. 

-135C.CCCC -44C.CCCC -C. -C. -C. -C. 

4 1 -2 920. CCCC -0. -C. -C. -C. -C. 

4 2 5080. CCCC 1940. CCCO -C. -0. -0. -C. 

4 3 -2350. COCO -330. COCO -C. -C. -C. -C. 

4 4 -74C.CCCC -1C. CCCC -C. -0. -C. -C. 

5 1 -385C.CCCG -0. -C. -C. -C. -C. 

5 2 -3840. CCCO -1430. CCCO -C. -C. -0. -C. 


5 

3 

-2340. CCCC 

760. CCCO 

-c. 

-c. 

-c. 

-c. 

5 


67C.CCCC 

190. CCCO 

-c • 

-c. 

-0. 

-c. 

5 

c; 

-IRC. CCCO 

ICO. CCCO 

-c. 

-c. 

-c. 

-c. 

6 

1 

lc20.CCCC 

-o. 

-c. 

-0. 

-c. 

-c. 

6 

2 

-3420. CCCO 

530. CCCO 

-c. 

-c. 

-c. 

-c. 

6 

A 

-123C.CC00 

-170. CCCO 

-c. 

-c. 

-c. 

-c. 

6 

4 

330. CCCC 

90. CCCO 

-c. 

-c. 

-c. 

-c. 

6 


330. CCCO 

320. CCCO 

-c. 

-c. 

-0. 

-0. 

6 

6 

5C.CCC0 

-40. CCCC 

-c. 

-c. 

-c. 

-c. 

7 

1 

-142C.CCCC 

-0. 

-C. 

-c. 

— C . 

-c. 

7 

2 

-54C.CCCC 

-19C.CCC0 

-c. 

-c. 

-c. 

-c. 

7 


2C.CCC0 

2060. OCGC 

-c. 

-c. 

-0. 

-c. 

7 

4 

2C5C.CCC0 

170. CCCO 

-c. 

-c. 

-0. 

-c. 

7 

5 

180. CCCO 

9C.CCC0 

-c. 

-C. 

-c. 

-c. 

7 

6 

-40. CCCC 

50. CCCO 

-c. 

-c. 

-0. 

-c. 

7 

7 

6C.CCC0 

20. CCCC 

-c. 

-c. 

-C. 

-c. 


J ones and Melotte (1953) Table V (Spheroidal earth) 


rv> ♦— * -C* 10 f\j vjj Nj 


Table VI (Cont) 


1 1 EPCCH= 1 96G . C 


JENSEN CAIN 


N N G 

2 1 3C411.2CC0 

2 2 2147. 4CCC 

3 1 2403. 5CC0 

-5125.3CC0 
-1338. 1CC0 
-3151. 8CC0 

6213. CCCO 
-2489. 8CC0 
-649.6CC0 
-4179. 4CC0 
-4529. 8CCO 
-2179. 5CCC 
5 4 7CC.8CC0 

5 5 -204.4CCC 

6 1 1625. 6CC0 

6 2 -344C.7CCC 

6 3 -1944. 7CC0 

6 4 -60.8C00 

6 5 277.5CC0 

6 6 69.7CC0 

7 1 -1952. 3CCC 

7 2 -485.3CC0 

7 3 321.2CCC 

7 4 2141. 3CC0 

7 5 105.1CCC 

7 6 22.70CC 

7 7 111.5CC0 


H 

GT 

HT 

GT T 

HTT 

-0. 

-C. 

-c. 

-c. 

-c. 

-5798. 9CC0 

-C. 

-c. 

-c. 

-c. 

-0. 

-0. 

-0. 

-c. 

-c. 

3312. 4CC0 

-C. 

-c. 

-c. 

-c. 

-157.9CC0 

-C. 

-c. 

-c. 

-c. 

-c. 

-C. 

-0. 

-c. 

-c. 

1487. CCGG 

-C. 

-c. 

-c. 

-c. 

-4C7.5CC0 

-c. 

-c. 

-c. 

-c. 

2 1 . C COO 

-c. 

-c. 

-c. 

-c. 

-0. 

-C. 

-0. 

-0. 

-c. 

-1182. 5CC0 

-c. 

-c. 

-0. 

-c. 

1CCC.6CC0 

-C. 

-c. 

-0. 

-c. 

43. CCCO 

-c. 

-c. 

-0. 

-c. 

138.5CC0 

-0. 

-c. 

-0. 

-c. 

-0. 

-G. 

-c. 

-0. 

-c. 

-79.6GC0 

-c . 

-c. 

-c. 

-c. 

-2C0.CCC0 

-c. 

-c. 

-0. 

-c. 

459.7CC0 

-c . 

-c. 

-c. 

-c. 

242. 1CG0 

-c. 

-c. 

-c. 

-c. 

-121.8C00 

-c . 

-c. 

-c. 

-c. 

-0. 

-c. 

-c. 

-c. 

-c. 

-575.8CC0 

-c . 

-0. 

-0. 

-0. 

-873.5CCC 

-c. 

-c. 

-c. 

-c. 

-34C.6CC0 

-c. 

-c. 

-c. 

-c. 

-11.8C00 

-c . 

-c. 

-0. 

-c. 

-111.6CC0 

-c. 

-0. 

-G. 

-c. 

-32.5CC0 

-0. 

-0. 

-0. 

-c. 


Jensen and Cain (1962) 


Table IV (Cont) 


DPI 74=P21*DP164+DP21*P1 64-0. 24026696*DP 154 611 

PI 75=P2 1 *P16 5—0 .2 324 8053* P 1 55 612 

DP 175=P21*DP 165+DP2i*P 165-0. 23248053*DP 155 613 

P 176=P2 1 *P 166-0. 22246941 *P 156 614 

DP176=P21*DP166+DP21*P 166-0. 2224694 1*DP 156 615 

PI 7 7=P2 1»P167— 0. 2 1023359*P157 616 

DP177=P21*DP167+DP2 i*P167— 0 .2 1023 359* DPI 57 617 

P178=P21*P 168-0. 19577308*P 158 618 

DPI 78=P2 I*DP168+DP21*P 168-0. 19 57730 8*DP 158 619 

P 179=P21*P 169-0. 17908787*P 159 620 

DP179=P21 *DP169+DP2 1*P 1 69— 0 . 17908787*DP159 621 

P1710=P21*P16 10-0. 160 17797*P 1510 622 

DP1710=P21*DP 1 610+DP2 l*P161 0-0. 160177 97*DP1510 623 

P1711=P21*P161l-0.13904338*P1511 624 

DP17ll=P21*DP1611+ DP 21*P161 1-0. 139043 38*DP1511 625 

P1712=P21*P1612-0.1 1 5 68409* PI 512 626 

DPI 712=P21*DP1612+DP2 l*P1612-0. 11568409*DP1512 627 

P1713=P21*P1613— 0.09010011*P1513 628 

DP1713=P21*DP1613+DP2l*PL613-0.09010011*DP15l3 629 

PI 7 14=P2 l *P16 1 4— 0 . 06 229143* P1514 630 

DPI 7 14=P21 *DP1614+DP2 l*P16l4-0 .06229143*DP1 514 631 

P1715=P21*P1615-0. 03225 80 6*P1515 632 

DP1715=P21*DP1615+DP21»P161 5- 0.032258 06*DP1515 633 

P1716=P21*P1616 634 

DP1716=P21 *DP 1616+DP21*P1616 635 

P171 7=P22*P1616 636 

DPI 717=16. 0*P1716 637 

AOR=AQR*AR 638 

C2=G( 17,2)*CP2+H( 17,2)*SP2 639 

C3=G( 17,3)*CP3+H(17,3)*SP3 640 

C4=G( 17,4)*CP4+H( 17,4)*SP4 641 

C5=G( 17»5)*CP5+H{ 17»5)*SP5 642 

C6=G( 17,6)*CP6+H{ 17,6)*SP6 643 

C7=G( 17,7 ) *C P7 + H ( 17» 7 )»SP7 644 

C8=G( 17,8) *CP8+Hl 17,8 )*SP8 645 

C9=G(17,9)*CP9+H{ 17,9)*SP9 646 

C10=G(17,10)*CP10+H(17,10)*SP10 647 

Cl 1=G( 17,11) *CP 1 1 +H( 17,11)*SP11 648 

C12=G ( 17 , 12) *CP12+H( 17, 12)* SP 12 649 

CL3=G( 1 7, 13)*CP13+H( 17, 13)*SP13 650 

C14=G< 17,14) *CP14+H( 17, 14)*SP14 651 

C 1 5=G ( 17,15) *CP 1 5+H ( 17,15)*SP15 652 

C16=G( 17,16) *CP16+H{ 17,16) *SP16 653 

C17=G( 17, 17)*CP17+H( 17, 17)*SP17 654 

BR=BR— 17 . 0*A0R* ( G ( 17 , 1 ) *P 1 7 1+C 2*P1 72+C3*P 173+C4*P 1 74+C5*P175+C6*P 1 655 

176+C7*P177+C8*P1 78+C9*Pl 79+C 10*P 17 10+C 11*P1711+C12*P1712+C13*P1713 656 

2 + C14*P1714+C15*P1715-*-C16*P1716+C17*P1717) 657 

BT = BT + AOR* ( G ( 17, 1 ) *DP171+C2*DP1 72+C3*DP 173+C4*0P 1 74+C5*DP175+C6*DP 658 

1176*C7*DP177+C8*DP178+C9*DP179+C10*DP17 10+C 11*DP1711+C12*DP1712+C1 659 

23*DP171 3+C 14*DP 17 14+C15*DP1715+C16*DP1716+C 1 7* DP 1717) 660 

BP=3P-A0R*( CGI 17,2) *SP2-H( 1 7 , 2 ) *CP2 ) *P 1 72 + 2. 0* ( G ( 1 7 , 3 ) * SP3-H ( 1 7 , 3 ) 661 

1*CP3) *P 173+3.0* (G( 17,4 >*SP4-H( 17,4 ) *CP4 ) *P1 74+4. 0* ( G ( 17 ,5)*SP5-H( 1 662 

27,5 ) *CP5 )* PI 75 + 5.0*1 G( 17,6) *SP6-H{ 17 , 6 ) *CP6 ) *P 1 76+6. 0* ( G ( 17 , 7 ) *SP7 663 

3-HC 17,7) *CP7)*P177+7.0*(G( 17 , 8 ) * SP8-H { 1 7 , 8 ) *CP8 ) *P1 78+8 . 0* ( G ( 17 , 9 ) 664 

4*SP9-H ( 17, 9) »CP9) *P1 79+9.0* (GC 17, 10 ) *SP10-H ( 17, 1 0 ) *CP 1 0 ) *P 1 7 10+ 10 . 665 

50* (G( 17, 11)*SP1 1 — H( 17,ll)*CPll)*P171l + 11.0*(G(17,12) * SP12— H (17,12) 666 

6*CP12)*P1712+12.0*(G( 17, 13) *SP13-H( 17, 13) *CP 13) *P 1713+13.0* (GC 17,1 667 

74) * S P 1 4— H ( 17, 14)*CP14)*P1714 + 14.0*(G( 1 7, 15 ) * SP1 5-H ( 1 7 , 1 5 ) *CP 1 5 ) *P l 668 

871 5+15. 0*(G( 17,16 )*SP 16-H( 1 7, 1 6) *CP16 ) *P 17 16+ 16. 0* ( G ( 17 , 1 7 ) * SP1 7-H 669 

9(17, 17)*CP17)*P1717) 670 

IF ( NMAX— 1 7 ) 18,18,17 671 



Table IV (Cont) 


C N= 18 672 

17 SP18=SP2*CP17+CP2*SP17 673 

CP 18=CP2*CP1 7— SP2*SP17 674 

P181=P21*P171-0.25024437*P161 675 

0P181=P2I*DP1 71+DP2 1*P 171-0. 25024437*DP 161 676 

P 1 82 = P2 1 *PI72— 0 .24926686* PI 62 677 

DP 182=P21*DP172+DP21*P 172-0. 24926686*DP 162 678 

P183=P2 L *P 1 73— 0.2463343 1*P 16 3 679 

DP183 = P2 l*DP173+DP21*P-173—0. 24633431* DP 163 680 

P184=P21*P174-0.24144672*P164 681 

DP184=P21*DP174+DP21*P1 74-0 . 24144672*DP 164 682 

P185=P21*P17 5-0.2 34604 10*P 165 68 3 

DP l 85=P2 1*DP1 75+OP2 1*P 175-0. 234604 10*DP 165 684 

P 1 86=P21 *P1 76-0. 22580645»P 166 685 

DP 186=P2i*DP176+DP21*P 176-0. 22580645*DP 166 686 

P 1 87=P2 1 *P177-0 .2 l 505376*P167 687 

DP187=P21*DP177+DP21*P177-0.21505376*DP167 688 

P188=P21*P178-0.20234604*P168 689 

DP188-P21*DP 1 78+DP2 1 *P 178— 0. 20234604*DP 168 690 

P 1 89=P21*P 179-0. 18768328*P 169 691 

DP189=P2 i*DP179+DP 21* Pi 79-0. 18768 328*DP 169 692 

P1810 = P2 l*P1710-0. 17106549*P1610 69 3 

DP1810=P21 *DP l 710+DP2 l*P1710-0. 17 1065 49* DP 1610 694 

P1811=P21*P1711—0. 15249 266* P1611 695 

0P1811 = P21*DP1711 + DP21»P1711-0. 152492 66*DP1611 696 

P1812=P21*P1712-0. 1 31 964 80* P 1612 697 

DP1812 = P21*DP1712 +DP2 l*P1712-0. 131 96480* DPI 6 12 698 

P1813=P21»P171 3-0 . 109 48191 *P 1613 699 

DP1813=P2l*DP1713*DP21*P1713-0. 109481 91 *OP 16 13 700 

P1814=P21*P1 714—0 .08504 398* PI 614 701 

OP 18 1 4=P21*DP1 714+DP2 1*P 171 4-0 .08504398*DP1 614 702 

P1815=P21*P1 715— 0 .05865 102*P 1615 703 

DP 1 8 1 5=P21*DP1 71 5+DP2 1*P 17 15-0 . 058651 02*0P161 5 704 

P1816=P21*P171 6— 0. 03030303*P 161 6 705 

DP1816=P21*DP1716+DP21*P1716-0.03030303*DP1616 706 

P1817=P21*P1717 707 

DP 1 8 1 7=P21 *DP171 7+DP2 1 *P 17 1 7 708 

P1813=P22*P1717 709 

DP1818=17.0*P1817 710 

AOR= AOR* AR 711 

C2 = G< 18,2)*CP2+H( 18,2)*SP2 712 

C3=G( 18,3)*CP3+H( 18,3)*SP3 713 

C4=G( 18,4)*CP4+H{ 18,4)*SP4 714 

C5=G( 18»5)*CP5+H(18»5)*SP5 715 

C6=G( 18,6)*CP6+H( 18,6)*SP6 716 

C7=G( 18, 7)*CP7+H{ 18,7)»SP7 717 

C8=G( 18, 8)*CP8+H( 18,8 ) *SP8 718 

C9=G( 18,9)*CP9+H( 18,9)*SP9 719 

C10=G( 18,10 )*CP10+H( 18,10 )*SP 10 720 

C11=G(18,11)*CP1 i+H( 18,11)*SP11 721 

C 1 2*G ( 18, 12 )*CP12+H( 18, 12)*SP12 722 

Cl 3=G( 18,13) *CP13+H( 18,13)*SP13 723 

C14=G ( 18,14) *CP 14+H ( 18,14)*SP14 724 

Cl 5 -G ( 18, 15)*CP15+H( 18, 15)*SP15 725 

C 1 6=G ( 18,16 )*CP16+Hl 18, 16)*SP16 726 

C 1 7 =G { 18,17) *CP 1 7+H { 18, 17)*SP17 727 

C 18 =G ( 18,18)*CP18+H(18, 18)*SP18 728 

BR = BR-18.0*A0R* ( G( 18,1 ) *P 18 1 + C 2*P1 82 +C3*P,183+C4*P1 84+C5*P 185+C6*P 1 729 

186+C 7*P187+C8*P188+C9*P189+C10*P1810+C11*P1811+C12*P1812+C13*P1813 730 

2+C14*P1814+C15*P1815+C16*P1816+Cl7*P1817+C18*P1818) 731 

BT=8T+AQR*CG( 18 , 1 ) *DP 181+C2*DP 182+C3*DP 1 83*C4*DP1 84*C5*DP 18 5+C6*DP 732 



Table IV (Cont) 


1186+C7*DP187+C8*DP188+C9*DP189+C10*DP1810+C1 1 *DP 1 81 l+C 12*DP 1812+C l 733 

23*DP1813+C14*DP1814+C15*DP1315+C16*DP1816+C17*0P1817+C18*DP1818) 734 

BP=BP— AOR* ( (G( 18, 2 )*SP2-H( 13,2 )*CP2) *P1 82+2.0* <G( 18,3) *SP3-H< 18,3) 735 

l*CP3)*P183+3. 0*10(18, A ) *SP4-H ( 1 8, 4 ) *CP4 ) *P1 84+4. 0* ( G ( 18,5)*SP5-H< 1 736 

28,5)*CP5)*P185+5.0*(G( 18 , 6 ) *SP 6-H< 1 8, 6 ) *CP6 ) *P 1 86 + 6. 0* 1 G ( 1 8 , 7 ) * SP7 737 

3-H(13,7)*CP7)*P187+7.0*(G(18,8)*SP8-H( 18, 8 ) *CP8)*P188+8.0*( G( 18,9) 738 

4*SP9— H( 18,9)*CP9)*P189+9.0*(G( 18 , 10 ) * SP 1 0-H ( 18 , 10 ) *CP10 ) *Pi 8 10 + 10 . 7 39 

50* (G( 18, 11 )*SPll-H( 18, 11 )*CPll ) *P1 81 1+11.0* (G( 18, 1 2 )* SP12-H ( 1 8 , 12) 740 

6*CP12)*Pl8l2+12.0*(Gl 18, 13)*SP13-H( 18 , 1 3 ) *CP 1 3 ) *P 18 1 3+1 3. 0* ( G ( 1 8, 1 741 

74)*SP14-H{ 18, 14 ) *CP14 ) *P 1 314+14.0* (G( 18, 15 )*SP15-H( 18,15 )*CP15)*P1 742 

8815+15. 0*( G( 18,16)*SP16-H( 18, 1 6) *CP16 ) *P 1816+16. 0* ( G ( 18 , 1 7 ) *SP1 7-H 743 

9( 18, 17)*CP17)«P1817+17.0*(G( 18, 1 8 ) *SP 1 8-HU8 , 1 8) *CP1 8 ) *P1 81 8 ) 744 

18 BP=BP/P 22 745 

B=SQRTIBR*BR+BT*BT+BP*BP) 746 

RETURN 747 

END 748 


Table V 


SUBROUTINE- CUNVRT (G,H,NMAX,K) 

K= I CONVERTS SCHMIDT TO GAUSS. K = 2 CONVERTS GAUSS TO SCHM 10 T 
DIMENSION G( I8*I8)«n{ 1 8 , 1 6 ) , S ( 16 « 18,2) 

S C 1 , 1 , 1 ) =- 1 . 

00 I N=2,NMAX 

S(N, l,l)=S(N-l,l, 1 )« FLOAT { 2* N-3)/ FLOAT (N-l) 

J = 2 
DU 1 

S ( N , M * I ) = S ( N » M— 1 » 1 )» SORT { FLOAT ( { N-M+ 1 )* J ) /FLOAT ( N+M-? ) ) 

J = I 

DO 2 N=2 t NMAX 
on 2 M=1 # N 

SIN, M, 2)=1 ./SlN.M, l) 

C.(N,MJ=G(N,M)*S(N,M f K) 

H(N,M)=H(N,M)*S(N,M,K> 

RE l URN 
END 


Table VI 


1 C EPCO* 1958.0 


ACAfS ET. A L . 


N 

2 

2 

3 

3 

3 

4 
4 
4 

4 

5 
5 
5 
5 

5 

6 
6 
6 
6 
6 
6 
7 
7 
7 
7 
7 
7 
7 


* 

G 


GT 

HT 

GTT 

HTT 

1 - 

3C6 IC • CCCC 

-c. 

11. CCCC 

-C. 

-c. 

-c. 

2 

-2C9C.CCC0 

5730. CCCO 

6 . CCCC 

l.CCCC 

- 0 . 

-c. 

1 

-1380. CCCC 

- 0 . 

-22. CCCC 

- 0 . 

- 0 . 

-c. 

2 

2910. CCCC 

-leec.ccco 

0 . 

-12. CCCO 

-c. 

-c. 

2 

189C.CCC0 

-4C.CCC0 

5. CCCC 

-2C.CCCC 

- 0 . 

-c. 

1 

95C.CCC0 

- 0 . 

3. CCCC 

-C. 

- 0 . 

- 0 . 

2 

-174C.CCCO 

-600. COCO 

-12.C0CC 

-l.CCOO 

- 0 . 

- 0 . 

3 

111C.CCCO 

370. COCO 

0 . 

5. CCCC 

- 0 . 

- 0 . 

4 

66C.CCC0 

30-C0C0 

-4. CCCC 

-e.cccc 

-c. 

-c. 

1 

113C.CCCC 

- 0 . 

-2. CCCC 

-c. 

-c. 

-c. 

2 

95C.CCCC 

260.CCQ0 

4. CCCC 

4.CC0C 

- 0 . 

- 0 . 

a 

620. CCCO 

-340. CCCO 

-2. CCCC 

-l.CCCC 

-c. 

-c. 

A 

-110. CCCC 

-40. CCCO 

l.CCCC 

l.CCCC 

-c. 

-c. 

c 

41C.CCC0 

-330. CCCO 

-2. CCCC 

-3. CCCC 

- 0 . 

-c. 

1 

-130.CCC0 

- 0 . 

-c. 

-c. 

- 0 . 

- 0 . 

2 

ICC. CCCC 

-380. COCO 

- 0 . 

-c. 

- 0 . 

-c. 

3 

17C.CCCC 

110. CCCO 

-c. 

-c. 

-c. 

-c. 

4 

-2CC * CCCC 

-130.CC00 

-c. 

-c. 

- 0 . 

-c. 

c 

-44C.CCC0 

-160. CCCO 

-c. 

-C. 

- 0 . 

- 0 . 

6 

-33C.CCC0 

40. CCCO 

-c. 

-G. 

-c. 

-c. 

1 

-12C.CCCC 

- 0 . 

-c. 

-c. 

-c. 

-c. 

2 

12C.CCC0 

40. CCCO 

— C • 

-c. 

- 0 . 

-c. 

3 

20. CCCC 

290.CC00 

-C. 

- 0 . 

- 0 . 

- 0 . 

4 

-380.CCC0 

30. CCCO 

-c. 

- 0 . 

-c. 

-c. 

c 

2C.CCC0 

2CC.C0C0 

-c . 

-c. 

- 0 . 

-c. 

6 

-25C.CCC0 

360. CCCO 

-c. 

- 0 . 

- 0 . 

-c. 

7 

-16C.CCC0 

70. COCO 

-C. 

-c. 

- 0 . 

-c. 


Adam , et al (1962) Table 4 

4dam et al (1963) Secular change from Appendix for 61 observatories 1954-1959 



Table VI (Cont) 


G 

C 

EPCCH- 

196G.C 

fiPRIL 64 

COEFFICIENTS 



N 

M 

G 

H 

GT 

h r 

GT I 

HTT 

2 

1- 

304 26 * 37 18 

0. 

18.9272 

0. 

-c. 

-c. 

2 

2 

-2173.5851 

5761.1544 

7.3484 

-1.8917 

-c. 

-c. 

3 

1 

-1547.8441 

0. 

-24.81C7 

C. 

-c. 

-c. 

3 

2 

2599.8304 

-1949.1511 

-0.8373 

-14.0435 

-0. 

-0. 

3 

3 

1574.2527 

201.4746 

0.7795 

-17.6824 

-G. 

-c. 

4 

1 

1322.6534 

0. 

-0.4255 

0. 

-c. 

-c. 

4 

2 

-2C08.9C58 

-441.7920 

-1C. 4983 

1.926C 

-c. 

-c. 

4 

a 

1275. C282 

232.8557 

3.3926 

3.9573 

-0. 

-c. 

4 

4 

877.2960 

-117.9688 

- 1.943C 

-8.9673 

-c • 

-c. 

5 

1 

956.7080 

0. 

0.8136 

C. 

-c. 

-c. 

5 

2 

796.9928 

149.2406 

5.3684 

-0.9192 

-c. 

-c. 

5 

3 

527.2557 

-266.3760 

-1.9078 

-1.6638 

-c. 

-c. 

5 

4 

-4C0.3606 

-3.6335 

-0. 156C 

3. 1923 

— c • 

-c. 

5 

c 

272.6586 

-262.0834 

C .8262 

-5.4794 

-c. 

-c. 

6 

1 

— 241.C056 

0. 

3.5132 

0. 

-c. 

-0. 

6 

2 

352.9445 

0.1744 

-0.6879 

1.8192 

-c. 

-c. 

6 

3 

23 1. 1537 

124. 3802 

2.4824 

2.9407 

-c. 

-c. 

6 

4 

-33. 1023 

-1C3.5910 

C.6225 

-0.7699 

-c. 

-c. 

6 

c. 

-146.6943 

-98.2729 

C.C353 

-0.3726 

-c. 

-c. 

6 

6 

-78.6225 

75.4376 

1.5714 

-0.2396 

-0. 

-c. 

7 

1 

57.6125 

-C. 

-c. 

-0. 

-G. 

-c. 

7 

2 

7G.6926 

5.9803 

-c. 

-c. 

-0. 

-c. 

7 

3 

20.4417 

85.5440 

-c. 

-c. 

-0. 

-c. 

7 

4 

-240.9729 

58.2674 

-c. 

-0. 

-c. 

-0. 

7 

c 

-19.1271 

-18.0492 

-c. 

-c. 

-c. 

-c. 

7 

6 

— C .4148 

-24.8694 

-c. 

-c. 

-c. 

-c. 

7 

7 

-1CC.4C9C 

4.5572 

-c. 

-c. 

-c. 

-c. 

8 

1 

89.9445 

-0. 

-c. 

-c. 

-G. 

-c. 

8 

2 

-47.2579 

-51.1593 

-c. 

-0. 

-0. 

-c. 

8 

'j 

-2. 1939 

-22.46CC 

-c. 

-c. 

-0. 

-c. 

8 

4 

-25.8317 

6.4347 

-c. 

-0. 

-0. 

-0. 

8 

*; 

-11.3775 

-37.7228 

-0. 

-0. 

-0. 

-c. 

8 

£ 

26. C820 

43.7049 

-c. 

-c. 

-0. 

-c. 

8 

7 

5.631C 

-3.6613 

-c. 

-c. 

-c. 

-c. 

8 

t 

5.8776 

-26.2983 

-c. 

-0. 

-0. 

-c. 


Daniels and Cain (1964) 
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Table VI (Con t) 


1 C E P C C H = 1955. C 


FINCH LEAICN 


N H C 

2 1-3C55C.CCCG 

2 2 “2270. CCCO 

3 1 -152C.CCCC 
303C.CCC0 
1580. CCCO 

nec.ccco 

“1510. CCCO 
1260. CCCO 
51C.CCC0 
55C.CCC0 
8CC.CCC0 
58C.CCC0 
-380.CCCC 
3 1C . CCCO 
-270. CCCO 
32C.CCC0 
2CC.CCC0 
“40. CCCO 
“150. CCCO 
-7C.CCCO 
ICC. CCCO 
50. CCCO 
20. CCCO 
-24C.CCC0 
-3C.CCCC 
0 . 

7 7 -HO. CCCO 


1 

GT 

HT 

GT T 

HTT 

-0. 

1C.CCCC 

“0. 

-C. 

“C. 

55C0.0CC0 

1C.CCCC 

“0. 

“0. 

“C. 

“0. 

-2C.5CCC 

“0. 

“0. 

“C. 

1900. COCO 

C. 

-16.5C00 

-0. 

“C. 

240. CCCO 

-C. 

-13.CCCC 

“0. 

“C. 

“0. 

-C.6667 

-c. 

“0. 

“C. 

-450. 00C0 

-12.6667 

7.3333 

“0. 

-0. 

290. CCCO 

1.3333 

4.CC0C 

“C. 

“C. 

-90. CCCO 

-2.2333 

-IC.CCCC 

“C. 

“C. 

“0. 

1.75CC 

“C. 

“0. 

-C. 

150.0CC0 

3.CCCC 

-2.5C00 

-0. 

“0. 

-310. CCCO 

-3.25CC 

-2.CC0C 

-0. 

“0. 

-40. CCCO 

l.CCCC 

1.5C00 

“0. 

“C. 

-170.0CC0 

1.25CC 

-6. CCCO 

“0. 

“C. 

-0. 

-1.60CC 

-C. 

“0. 

-0. 

20. COCO 

C.2CCC 

5.2CCC 

“C. 

“C. 

ICO. COCO 

2.80CC 

1.2CCC 

“0. 

“C. 

—5C.0CC0 

0.2CCC 

“2.6C0C 

“0. 

“0. 

- 140. OCCO 

“l.COCO 

1 • 6000 

“0. 

“0. 

90. COCO 

0.4CCC 

0.4CCC 

-c. 

-c. 

“0. 

C • 1667 

-0. 

“0. 

“C. 

-20. OCCO 

l.CCCC 

-3.3333 

“0. 

“C. 

110. OCCO 

1.1667 

-0.8333 

-0. 

-c. 

0. 

l.CCCC 

2.CC0C 

“0. 

“C. 

“1C.OCOO 

-C.5CCC 

-1.667C 

-0. 

“C. 

“30. OCCO 

0.1667 

2.5CC0 

-0. 

“C. 

-10.0000 

-1.6670 

-0. 

“C. 

“C. 


Finch and Leaton (I960) 
Lea ton (1962) 



Table VI (Cont) 


I C E PCCH= 1960.0 FOLGERE 


N 

N 

G 

h 

GT 

HT 

GT T 

HIT 

2 

1- 

305C9.CCC0 

0. 

-c. 

-c. 

-c. 

-c. 

2 

2 

-2181. CCCO 

5841. COCO 

-c. 

-c. 

-C. 

-c. 

3 

1 

-1464. CGCO 

0. 

-c. 

-C. 

-0. 

-c. 

3 

2 

2971. CCCO 

-1988. CCCO 

-0. 

-c. 

-c. 

-c. 

3 

3 

1673. CCCO 

198. CCCO 

-c. 

-c. 

-0. 

-c. 

4 

1 

1147. CCCO 

0. 

-c. 

-0. 

-0. 

-c. 

4 

2 

-2012. CCCO 

-390.0CC0 

-0. 

-C. 

-C. 

-c. 

4 

a 

1151. CCCO 

291. CCCO 

-c. 

-c. 

-c. 

-c. 

4 

4 

831. CCCO 

-174. CCCO 

-C. 

-c. 

-c. 

-c. 

5 

1 

997.CCCG 

0. 

-0. 

-0. 

-0. 

-c. 

5 

2 

863. CCCO 

173. CCCO 

-0. 

-0. 

-c. 

-C. 

5 

TJ 

583. CCCO 

-292. CCCO 

-C. 

-c. 

-c. 

-c. 

5 

4 

-312. CCCC 

69. CCCO 

-c. 

-c. 

-c. 

-c. 

5 

c 

237. CCCO 

-206. COCO 

— c . 

-0. 

-0. 

-c. 

6 

1 

-267. CCCO 

0. 

-0. 

-o. 

-0. 

-c. 

6 

2 

2 72.CCCC 

-39. COCO 

-c. 

-0. 

-c. 

-c. 

6 


251.CCCC 

161. CCCO 

-c. 

-c. 

-C. 

-c. 

6 

4 

-52. CCCO 

-127.0CC0 

-c. 

-c. 

-0. 

-C. 

6 

t. 

-1 13.CCCC 

-147. COCO 

-c. 

-0. 

-c. 

-c. 

6 

6 

-8. CCCO 

88. CCCO 

-c. 

-C. 

-c. 

-c. 

7 

1 

48. CCCO 

0. 

-c . 

-c. 

-c. 

-c. 

7 

2 

131 .CCCC 

1 .CCOO 

-c. 

-C. 

-0. 

-0. 

7 

a 

-15. CCCO 

62. CCCO 

-c. 

-c. 

-c. 

-c. 

7 

4 

-224. CCCO 

38. CCCO 

-c. 

-c. 

-c. 

-c . 

7 

c 

13. CCCO 

19. CCCO 

-c. 

-0. 

-c. 

-c. 

7 

6 

-71. CCCC 

4. COCO 

-c. 

-C. 

-C. 

-c. 

7 

7 

-67. CCCO 

-57. CCCO 

-c. 

-c. 

-c. 

-c . 

8 

1 

1C4.CCC0 

0. 

-c . 

-c. 

-c. 

-c. 

8 

2 

-64. CCCO 

-69 .CCCC 

-c. 

-c. 

-0. 

-c. 

8 

a 

8. CCCC 

7. CCCO 

-c. 

-C. 

-0. 

-C. 

8 

4 

-18. CCCC 

37. CCCO 

-c. 

-0. 

-C. 

-c . 

8 


-42. CCCO 

22. COCO 

-c. 

-c. 

-C. 

-c. 

8 


55. CCCC 

8. CCCO 

-c. 

-c. 

-C. 

-c. 

8 

7 

-6. CCCO 

62.0CC0 

-0. 

-C. 

-0. 

-c. 

8 

€ 

-33. CCCO 

3. CCCO 

-c. 

-c. 

-c. 

-c. 

9 

1 

-16. CCCO 

C. 

-c . 

-c. 

-C. 

-c. 

9 

2 

6. CCCC 

-6. CCCO 

-c. 

-c. 

-0. 

-c. 

9 

a 

13. CCCO 

-21. COCO 

-c. 

-c. 

-0. 

-c. 

9 

4 

-5. CCCO 

-l.CCCO 

-C. 

-c. 

-0. 

-c. 

9 

5 

3 3 . C C C C 

-70. CCCO 

-c . 

-c. 

-C. 

-c. 

9 

6 

-1C. CCCC 

5. CCCC 

-c. 

-0. 

-0. 

-c. 

9 

7 

-7. CCCC 

-ll.OCCO 

-c. 

-c. 

-c. 

-c. 

9 

£ 

61. CCCC 

0. 

-c. 

-c. 

-c. 

-c . 

9 

9 

C. 

0. 

-c . 

-c. 

-C. 

-c. 


Fougere (1964) 
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